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ABSTRACT 

The optical afterglow of long-duration GRB 071003 is among the brightest yet to be detected from 
any GRB, with R w 12 mag in KAIT observations starting 42 s after the GRB trigger, including 
filtered detections during prompt emission. However, our high S/N ratio afterglow spectrum displays 
only extremely weak absorption lines at what we argue is the host rcdshift of z — 1.60435 — in contrast 
to the three other, much stronger Mg II absorption systems observed at lower redshifts. Together with 
Keck adaptive optics observations which fail to reveal a host galaxy coincident with the burst position, 
our observations suggest a halo progenitor and offer a cautionary tale about the use of Mg II for GRB 
rcdshift determination. We present early through late-time observations spanning the electromagnetic 
spectrum, constrain the connection between the prompt emission and early variations in the light curve 
(we observe no correlation), and discuss possible origins for an unusual, marked rebrightcning that 
occurs a few hours after the burst: likely either a late-time refreshed shock or a wide-angle secondary 
jet. Analysis of the late-time afterglow is most consistent with a wind environment, suggesting a 
massive star progenitor. Together with GRB 070125, this may indicate that a small but significant 
portion of star formation in the early universe occurred far outside what we consider a normal galactic 
disk. 

Subject headings: gamma rays: bursts — gamma-ray bursts: individual: 071003 



1. INTRODUCTION 

Concurrent observations of long- wavelength afterglow 
and ongoing gamma-ray burst (GRB) activity should, in 
principle, yield important constraints on th e nature of 
the p hysical processes of the emission (eg., iKobavashl 
2000). However, as a GRB typically lasts less than 
100 s, it is challenging for large ground-based opti- 
cal/infrared follow-up facilities to react to a GRB alert 
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quickly and take data during the prompt phase. Multi- 
color observations, which provide vital information on 
the emission mechanism, are even more difficult to ob- 
tain during the prompt phase because of the added 
overhead associated with changing filters. Nevertheless, 
due to the coordinated effort s of r ecent space missions 
(HETE-II. Hicker et al.ll200l : Swift. lOehrels et al.ll200l 
to detect GRBs and various ground-based optical follow- 
up programs, observations during the prompt phase of 
GRBs are no longer uncommon — the optical after- 
glows (OAsj_of_severaIdozen_Cj^^ 
(og JAkerlof et al.ll999tlVestrand et al.l2006HYost et all 
2003) during gamma-ray emission, and multi-color opti- 
cal data have been obtained in a handful of cases (e.g., 
iBlake et al.ll2005l: iNvsewander et al.ll2007f) . 

Observations of GRBs in the past several years 
have also revealed a rich demography in OA behav- 
ior. Some OAs have monotonic power-la w decays (e.g., 
iLi et alj l2003bl : lLaursen fc Stanekl [2 003). while others 
have plat eau (e.g., iRvkqff et al] I2006D and rebrighten- 
ing (e.g., Wozniak ct al. 2006) phases. Even among 
GRBs with relatively simple behavior, however, short- 
timescale features not predicted in the basic shock mod- 
els often appear in sufficiently well-sampled data. Var- 
ious modifications to the standard picture have been 
proposed to explain such observatio ns, including the 
presence of a jet with single (e.g., ISari et aLI Il999f ) 
or multiple fe. g.. iBerger et ail I2003D components, re- 
fres hed shocks (|Zhang e t al. 200 6|), central engin e activ- 
ity (jKocevski et al.l 120071; IChincarini et alJ[2007h . gravi- 
tational microlensing (jGarnavich et al.ll2000D. and den- 
sity i rregularity in the GRB environment ( Hollan d et al] 
2003). Observationally, constraints on the change in 
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the afterglow color and the spectral energy distribution 
(SED) play an important role in limiting the viability of 
models for a particular GRB. 

The question of the nature of the GRB itself is in- 
timately tied to the question of its environment and 
origins. At interme diate to late times, spectroscopy of 
the a fterglow (e.g., iProchaska et all 12001 iD'Elia et all 
20071) and deep imaging of the host en vironment (e.g., 
Bloom et ail 120021 iFruchter et~al]|2006ft can help estab- 
lish the nature of the GRB's progenitor and environment, 
connecting what we learn about the burst itself to the 
larger question of its origins and place in the early uni- 
verse. 

In this paper, we report on our photometric and spec- 
troscopic observations of GRB 071003 with various tele- 
scopes from the prompt phase to late times. In §2 we 
describe the observations, and in §3 we present the reduc- 
tions. The analysis of the light curves and the constraints 
on the changes in the colors and SEDs are given in §4. 
The conclusions, including the implications of the ex- 
tremely unusual spectrum of this event, are discussed in 
§5. We assume H Q = 71 km s _1 Mpc -1 , f_M = 0.3, and 
Oa = 0.7 throughout the paper. These reports should 
be considered the final analysis of our group's data on 
GRB 071003, superseding any previously announced re- 
sults (e.g., in the GCN Circulars). 




t-t,rigge, (Sec) 

Fig. 1.— Light curve from the Swift BAT of GRB 071003, with 
optical photometry from KAIT and P60, and the optical light curve 
model discussed in Section 4.3, overplotted. The GRB is dominated 
by a complicated, spiky emission episode in the first 30 s, but a 
pulse is also observed much later, at 150 s. Optical data points 
(all from KAIT, except one _R-band measurement from the P60), 
by contrast, show a power-law decay at early times followed by a 
slow-rising "bump." Here the V and I filtered observations have 
been offset to match the R and unfiltered points based on the 
relative colors at 2000 s. 



2. OBSERVATIONS 
2.1. BAT/XRT Observations 

On 2007 October 3, 07:40:55 UT (defined as t = in 
this paper; UT dates are used throughout), a bright GRB 
triggered the Burst Alert Tele scope (BAT) onboa rd the 
Swift satellite (trigger 292934 iSchadv et all [20071) . The 
first GCN notice was distributed within 16 s. Unfortu- 
nately, Swift was still returning to normal observations 
after its 2007 August gyro failure, but it did slew to the 
position after 22 ks and began observations using the 
X-Ray Telescope (XRT). 

We downloaded the Swift BAT and XRT data from 
the Swift Archive 17 and quicklook data site. 18 The XRT 
and BAT spectra were fitted using ISIS 19 . 

The XRT data were processed with version 0.11.4 
of the xrtpipeline reduction script from the HEA- 
soft 6.3. 1 20 software release. We employ the latest 
(2007 December 4) XRT calibration files. Our reduc- 
tion of XRT data from cleaned event lists output by 
xrtpipeline to science -ready light curves and sp ectra 
is described in detail bv lButler fc Kocevsk i (2007aj). We 
use the latest calibration files from the 2007 September 
24 BAT database release. We establish the energy scale 
and mask weighting for the BAT event mode data by run- 
ning the bateconvert and batmaskwtevt tasks. Spec- 
tra and light curves are extracted with the batbinevt 
task, and response matrices are produced by running 
batdrmgen. To produce the BAT spectra, we apply the 
systematic error corrections to the low-energy BAT spec- 
tral data as suggested by the BAT Digest Web site 21 , 
and fit the data in the 15-150 keV band. The spectral 

17 ftp:/ /legacy.gsfc. nasa.gov/swift/data 

18 \protect 'ittp: / /swift. gsfc.nasa.gov/cgi-bin/sdc/ql 



normalizations are corrected for satellite slews using the 
batupdatephakw task. 

The burst exhibits one dominant emission episode of 
duration dt « 30 s, followed by a minor pulse ~ 150 s 
later of duration ~ 20 s. The total duration is Tg = 
148 ± 1 s, 22 placing it clearly into the long GRB class. 
The primary pulse is resolved into multiple pulses. The 
gamma-ray light curve is shown in Figure [TJ overplotted 
with early-time photometry from KAIT and P60 (dis- 
cussed in §3.2 and §3.5, respectively). 

The time-integrated BAT spectrum from t = —10.3 to 
t = 169 s is acceptably fitted (x 2 /^ = 47.64/55, where 
v is the number of degrees of freedom) by a power-law 
model, with photon index a — —1.3 ± 0.1 and energy 



fluence 5 7 = (1.7 ±0.1) x 10 



" 5 erg cm" 2 (15-350 keV). 



The main emission episode (t = —1.4 s to t = 22.8 s) is 
harder {a = -1.08 ± 0.03, S 7 = (1.51 ± 0.03) x 10~ 5 erg 
cm~ 2 , x 2 l v — 56.71/55), while the final pulse (t = 131— 
169 s) is softer (a = -1.8 ± 0.2, 5 7 = 1.2±°, 1 x 10~ 6 erg 

cm" 2 , X 2 /v = 41.15/55). 

X-ray observations with the XRT began 6.2 hr af- 
ter the BAT trigger. The X-ray light curve measured 



until t 



5 x 10° s is well fitted by a power-law 



19 \protect ittp://space. mit.edu/CXC/ISIS 

20 \protect ittp://heasarc. gsfc.nasa.gov/docs/software/lheasoft/ 

21 \protect 'ittp: / /swift. gsfc.nasa.gov/docs/swift/analysis/bat_digest. 



time decay ^-i-ssio- 03 . The time-integrated spectrum 
is well fitted (x 2 / l/ — 48.47/54) by an absorbed power- 
law model [photon index T — 2.14 ± 0.12, unabsorbed 
F x = (5.8 ± 0.4) x 10~ 13 erg cm" 2 s" 1 ]. The equiva- 
lent H column density, 7V H = (2.2 ± 0.4) x 10 21 cm -2 , 
is marginally consistent with the expected Galactic col- 
umn d ensity in the source dire ction, TVjj = 1.1 x 10 21 
cm -2 (jDickev fc Lockmanlll990fl. Examining the X-ray 
hardness ratio fe.g.. iButler fc Kocevskilr2007bf ). there is 
no evidence for spectral evolution during the XRT obser- 
vation. 

22 All uncertainties quoted in this paper are la, except where 
hpamified otherwise. 
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2.2. KAIT Observations 

The Katzman Automatic Imaging Telescope (KAIT) 
is a 0.76-m robotic telescope at Lick Observatory that 
is dedicated to searching for and observing supernovae 
and monitoring other variable or ephemeral celestial 
phenomena. It is equipped with a Finger Lakes In- 
strument (FLI) ProLine PL77 back-illuminated CCD 
camera having a resolution of 0'.'8 pixel -1 and a to- 
tal field of view (FOV) of ~ 6.'8 x 6'.8. More in - 
formation on KAIT can be found in iLi et alJ (|2000f ), 
iFilippenko eTafl (|200ll ). and iFilippenkol pOOa) . while 
the KAIT GRB alert system is described in detail by 
ILi et all (|2003bh . Nota ble KAIT observ ations of GRBs 
include GRB 021 211 (|Li et all l2003eJ). GRB 051111 
(jButler et al.ll200l . GRB 060210 (IL1I21M Li et al. 200 8 
in preparation), and GRB 080319B (|Bloom et al.|[200l . 

Several improvements have been implemented for the 
KAIT GRB ale rt system since the description given by 
ILi et all (|2003b[ ). An FLI PL77 camera has replaced the 
Apogee AP7 camera, offering a much faster readout time 
(1.2 s for FLI vs. 11.0 s for Apogee). A new feature has 
been incorporated into the software so the system can 
easily terminate an ongoing exposure in preparation for 
the GRB response sequence. Most importantly, a real- 
time image-processing pipeline has been developed to 
compare the KAIT images to archival Digital Sky Survey 
(DSS) images to identify new objects. Astrometry solu- 
tions are derived for the KAIT images by matching the 
detec ted objects to the USNO Bl catalog (jMonet et al.1 
2003), providing coordinates to any new objects to a pre- 
cision of ~0".2. Point-spread-function (PSF) fitting pho- 
tometry is also performed on new objects, and calibrated 
to the red magnitudes of the stars in the USNO Bl cat- 
alog. The image-processing results are displayed in real 
time on a website. 23 

For GRB 071003, the KAIT GRB alert program re- 
ceived the GCN socket notice at t = 16 s. The system 
immediately terminated the ongoing supernova search 
program and began to slew the telescope to the GRB po- 
sition. After slewing from close to meridian to an hour 
angle of 4.2 hr, a sequence of 5x5 s unfiltered images 
began at t — 42 s. KAIT then switched to a sequence 
that alternated with 20 s V, I, and unfiltered images. 
Finally, the sequence converted to 20 s / and unfiltered 
images. Because of the physical west hour angle limit 
of 4.7 hr, KAIT only finished part of this pre-arranged 
sequence. In total, 56 images were obtained in the V, I, 
and unfiltered passbands from t = 42 to 1628 s, with full 
width at half-maximum intensity (FWHM) of ~3". 

Visual inspection of the image-processing results re- 
vea led a true new object, first reported by our group 
(|Lill2007D . measured at 12.8 mag at a position of a = 
20 /l 07 m 24 s .12, 5 = +10°56'51".8 (equinox 2000.0; ap- 
proximate lcr astrometric uncertainty 0".3). Our can- 
didate OA was subsequently confirmed by observations 
from the automated Palo mar 60-inch (1.5 m) telescope 
(P60: ICenko fe Foxll2007| ). Further preliminary analysis 
of OA early-time behavio r from the KAIT observations 
was reported by ILi et all (|2007). Figure [2] shows a se- 
quence of the KAIT images for the OA of GRB 071003. 
An 80" x 80" section is shown for the first and fifth unfil- 

23 \protect ittp://hercules. berkeley.edu/grbdata/grbfinder.gif . 



tered 5 s image and a 20 s unfiltered image that started 
at t = 431 s. As s een in Figure |2] and r eported by several 
groups (0120071 : ICenko fc Foi^OOTI iMisra et al.l 120071: 



ILi et al.ll2007l ). a bright (R ~ 11 mag) foreground star 
is located 6".5 west of the OA of GRB 071003. As dis- 
cussed in §3, the presence of this bright star complicates 
the photometry for the OA, and various methods have 
been used to minimize its contamination. 

2.3. P60 Observations 



Th e Palomar 60-inch telescope (P60; ICenko et al.l 
l2006h automatically responded to the Swift trigger for 
GRB 071003, beginning a pre-programmed sequence of 
observations at 07:43:51 UT (176 s after the trigger). 
Observations were taken in the Kron R, Sloan i' and z' , 
and Gunn g filters at large airmass (> 2.5). Individual 
images were reduced in real time by our automated re- 
duction pipe line. Manual inspe ction revealed a fading 
point source (|Cenko fc F ox 2007) in all four filters at the 
location reported byO (|2007h . 

A second epoch of observations was manually sched- 
uled for the night of UT October 4. In an attempt to 
lessen the contamination of the nearby bright saturated 
star, these observations were taken in the Johnson V- 
band filter in relatively short (30 s) exposures. A se- 
quence of 30 images was obtained. 

2.4. AEOS Observations 

The 3.6-m US Air Force Advanced Electro-Optical 
System (AEOS) telescope, located at the Maui Space 
Surveillance System on Haleakala 24 , observed the 
OA of GRB 071003 with the A EOS Burst Camera 
(ABC, iFlewelling-Swan et all 120061) . ABC has a back- 
illuminated 2048 x 2048 pixel EEV chip, with a scale 
of 0".189 pixel" 1 and a FOV of - 6f5 x 6.'5. Because 
there is no direct internet access to AEOS, after Swift 
detected the GRB, a FAX alert was automatically sent 
to the AEOS control room, to initiate a series of Target- 
of-Opportunity (ToO) observations. 

The AEOS observations of GRB 071003 are all unfil- 
tered 10 s exposures. The first batch of images started at 
~9 minutes after the BAT trigger, and 238 images were 
observed until t « 83 minutes, all with very good image 
quality (FWHM w 0".9). The second batch of images 
started at t « 205 minutes, and 56 images were observed 
until t k, 222 minutes. Due to the large airmass for these 
observations and the degraded seeing conditions, how- 
ever, the images have rather poor quality. We have tried 
various methods to measure the brightness of the OA in 
these images but failed. Accordingly, only the first batch 
of 238 images is analyzed in this study. Prelim inary anal- 
ysis of the AEOS observations is reported bv lSwan et al.l 
(l200l . 

2.5. Keck I/Gemini-S observations 

In response to the detection of the OA of GRB 071003, 
we organized a campaign to obtain spectroscopy and 
late-time photometry with the 10-m Keck I and the 8-m 
Gemini-S telescopes. At t ~ 2.6 hr, we attempted to 
observe the OA with the HIRES spectrograph at Keck 

24 Based on data from the Maui Space Surveillance System, 
which is operated by Detachment 15 of the U.S. Air Force Research 
Laboratory's Directed Energy Directorate. 
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t = 42 - 47 s 
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t = 67 - 72 s 




t = 431 - 451 s 




Fig. 2. — Sequence of KAIT images for the OA of GRB 071003. An 80" X 80" section is shown for the first and fifth unfiltered 5 s images 
and for a 20 s unfiltered image that started at t = 431 s. The OA is the central object in the circles. It is well detected in the early images 
and rapidly fades. The image quality is poor owing to the very high airmass of the object. 



I, but the data are of poor signal- to- noise ratio (S/N ) 
and no obvious lines were detected (|Perlev et al.ll2007bh . 
Just before the HIRES spectroscopy started, we also ob- 
tained guider images for the OA, providing important 
photometric coverage during a gap in the photometry 
obtained elsewhere (see §3.4). The guider images have a 
scale of 0".37 pixel" 1 with a FOV of 53". 5 x 71".3. 

On 2007 October 4, we observed the GRB 071003 OA 
with the Low R esolution Imaging Spectrometer (LRIS; 
lOke et all (|1995[ )1 on Keck I. Anticipating significant fad- 
ing of the OA, a series of deep 300 s images was taken 
with the g and R filters under excellent seeing conditions 
(FWHM k, 0".5). Inspection of the images reveals that 
the OA was still bright and saturated in most of the im- 
ages. Consequently, only a single image in each of the g 
and R bands, where the OA is not saturated, is analyzed 
in this study. LRIS uses a beamsplitter to separate the 
light between two arms, red and blue. Both the blue and 
red cameras have a usable FOV of ~ 6.'0 x 7.'8. The red 
camera used a back-illuminated Tek 2048 x 2048 pixel 
chip with a scale of 0".215 pixel -1 , while the blue cam- 
era has a mosaic of two 2048 x 4096 pixel Marconi chips 
with a scale of 0".135 pixel -1 . 

Encouraged by the brightness of the OA, we also per- 
formed LRIS spectroscopy of the OA. A preliminary 
analysi s of the spectrum is reported by iPerlev et al.1 
(|2007al) . and a more detailed analysis is presented in §3.5. 

We performed more LRIS imaging for the OA of 
GRB 071003 on 2007 October 8, 9, 10, 11, and 15, us- 
ing various combinations of u, g, V, and R filters. The 
presence of the very bright star presents a significant 
challenge to extracting useful data on the OA, as its 
diffraction spikes change positions and intensity accord- 
ing to the time and seeing conditions of the observations. 
Unfortunately, observations on 2007 October 8 were ad- 
versely affected by diffraction spikes and poor seeing, and 
were not usable. The data taken on 2007 October 15 are 
seriously affected by clouds, and do not provide an inter- 
esting limit to the brightness of the OA, so they are not 
used in this study. 

We also triggered our TOO program (GS-2007B-Q- 
2; PI H.-W. Chen) for GRBs with the Gemini-S tele- 
scope and obtained g-, r-, i-, and z-band images with the 
CMOS camera on 2007 Oct. 5 and 6. The CMOS camera 
is equipped with three back- illuminated EEV 2048 x 4608 
pixel chips. For our observations, the camera is used in 
a 2x2 binning mode with a scale of 0".146 pixel -1 and 
a FOV of ~5'.5x5'.5. Unfortunately, the 2007 October 



5 images are badly affected by bleeding from the very 
bright star and are not used in this study. 

As part of the efforts to follow the evolution of the 
OA of GRB 071003, we also performed adaptive op- 
tics (AO) observati ons with Keck I on 2007 October 19 
(jPollack et al.ll2007t) . The details of the AO observations 
can be found in §3.6. 

2.6. Radio Observations 

GRB 071003 was observed with the Very Large Array 
(VLA) 25 on various occasions. We made the observa- 
tion in the B configuration array. We used VLA source 
1950+081 as phase calibrator for 4.86 GHz (C) band ob- 
servations and 2001+104 for 8.46 GHz (X) band obser- 
vations. The data were analyzed using standard data re- 
duction routines of the Astronomical Image Processing 
System (AIPS). The first observation took place on 2007 
October 5 in the X band with flux density of 393 ± 55 
/zJy. Since then we made six observations in the X band 
and three observations in the C band (Table [6]). 

3. DATA REDUCTION 

The bright star in the neighborhood of the OA of GRB 
071003 makes it a challenge to measure reliable photom- 
etry from the data described in §2. In this section we 
describe the methods used to minimize its contamina- 
tion. 

3.1. Photometric Calibrations 

For photometric calibrations, the field of GRB 071003 
was observed in B, V, R, and / on two photometric 
nights (2007 October 7 and 8) at Lick Observatory, us- 
ing both KAIT and the Lick Nickel 1-m telescope. About 
a dozen Landolt standard-star fields (Landolt 1992) were 
observed at different airmasses throughout each photo- 
metric night. Photometric solutions to the Landolt stan- 
dard stars yield a scatter of ~0.02 mag for all the filters. 
The GRB 071003 field was also observed for several sets 
of BVRI images with different depth on both nights. 
The photometric solutions are used to calibrate a set of 
local standard stars in the GRB 071003 field. Because 
the GRB 071003 field is quite crowded, the number of cal- 
ibrated local standard stars is large (Table [T]). A finder 
for a subset of 23 relatively bright local standard stars 

25 The NRAO is a facility of the National Science Foundation, 
operated under cooperative agreement by Associated Universities, 
Inc. 
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Fig. 3. — Finder chart for a subset of local standard stars in 
the field of GRB 07f003. The field of view is 6'. 6x6'. 6. North is 
up and east is to the left. The displayed image is the unfiltered 
template taken with KAIT on 2007 October 9. 

is in Figure [3J As seen in Table [H the local standard 
stars in the field of GRB 071003 are well calibrated, with 
standard deviation of the mean (SDOM) of ~0.01 mag 
for all the BVRI bands. We refer to this calibration as 
the "Lick calibration" throughout the rest of the paper. 

Several Landolt standard-star fields were also observed 
with LRIS at Keck I: in the u, g, and R bands on 2007 
October 9, and in the V band on 2007 October 11. As the 
number of the observed standard-star fields is small, it 
is not possible to derive a complete photometric solution 
for either night. Since the GRB field was observed at 
similar airmasses with some of the standard-star fields, 
we can treat the LRIS filters as standard and derive the 
magnitudes for the local standard stars via differential 
photometry. Unfortunately, this procedure suggests that 
the 2007 Oct. 9 night was not photometric, as different 
standard-star observations yield somewhat different zero 
points. The 2007 October 11 night was photometric, but 
only the T^-band standard stars were observed. 

We elected to use the Lick calibration as the founda- 
tion for all the photometric calibrations, except in the 
case of the u band. The Lick-calibrated magnitudes are 
in BVRI, and can be reliably converted to t he q, r, and i 
band s using color transformation equatio ns dJester et al. 
2005). The conversion to the z band (Rodgcrs et al. 
20061) is somewhat problematic, and as a result we adopt 
a relatively large uncertainty for the converted magni- 
tudes. For the u band, only two standard-star fields 
were observed with LRIS on 2007 Oct. 9, and they give 
a difference of 0.30 mag in the zero points. We chose to 
calibrate the GRB 071003 field with the standard-star 
field that is closer in time of GRB observation, but we 
added an uncertainty of 0.30 mag to all the calibrated 
magnitudes. We note that the true error for the it-band 
calibration may be higher than 0.30 mag due to the non- 
photometric conditions on 2007 October 9. 

3.2. KAIT Data Reduction 
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Fig. 4. — Illustration of using image subtraction to remove the 
contamination of the bright nearby star to the OA of GRB 071003. 
The KAIT image subtraction code is demonstrated here, (a) An 
80" X 80" section of the original 20 s unfiltered KAIT image of the 
OA taken at t = 43f s; (b) the same section after image subtraction 
of the central 50" X 50" using an unfiltered template image after the 
OA has faded; (c) an 80" X 80" section of the combined unfiltered 
AEOS image at t = 5002.6 s; and (d) the same section after image 
subtraction of the central 30" X 30" using a hand-made template 
image. See text for more details. 

The KAIT data were automatically processed with bias 
and dark current subtraction and flat-fielding. The PSF 
of the OA is seriously affected by the bright star which 
is less than 10 pixels away in the KAIT images. Con- 
sequently, normal PSF-fitting photometry cannot fit the 
peak and background of the OA simultaneously to pro- 
duce a reliable measurement. 

We use image subtraction to remove the contamina- 
tion of the bright star. To generate template images for 
subtraction, KAIT imaged the GRB 071003 field in the 
unfiltered mode and in the V and / filters for the next 
several nights after the burst. To make sure the bright 
star is not saturated, short (5 s) exposures were used, 
and 50-100 images for each filter were acquired to en- 
sure high S/N in the combined images. As discussed in 
§4, the GRB OA was still reasonably bright in the sec- 
ond night after the burst, so we used the images obtained 
at 4-6 days after the burst as the template for the field 
without significant OA contribution. O ur image subtrac- 
tion c ode is based on the ISIS package (|Alard &; Luptonl 
1998) as modified by B. Schmidt for th e High- z Super- 
nova Search Team ( Schmid t et al.l[T998f ). An illustration 
of the image subtraction is presented in the top panels 
of Figure [2 

The Lick calibration was used to transform the KAIT 
instrumental magnitudes to the standard Johnson V and 
Cousins / passbands, with proper color terms measured 
from the photometric nights. We also find that the com- 
bination of the KAIT optics and the quantum efficiency 
of the FLI CCD camera makes the KAIT unfiltered ob- 
servations mostly mimic the R band. During the two 
photometric nights, unfiltered observations of the Lan- 
dolt standard-star fields were also performed. Analysis 
of these images indicates that the KAIT unfiltered mag- 
nitudes can be effectively transformed to the R band, 
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Fig. 5. — Early-time light curve of the optical afterglow of 
GRB 071003 using KAIT photometry, supplemented by observa- 
tions from P60 and AEOS. The gamma-ray light curve from the 
BAT is overplotted in gray (scaled arbitrarily). A clearly additive 
"bump" at 100-500 s is apparent. Photometric follow-up observa- 
tions continued after 2000 s with P60 and AEOS, as well as with 
Gemini and Keck in subsequent ni ght s; the complete 16-day optical 
light curve is presented in Figure HOI 

with a relatively large color term and an rms of ^0.05 
mag, s i milar to the earlier results we reported (|Li et alJ 
I2003bllal ; iButler et al.ll2006l ). 

To increase the S/N, the late-time KAIT images of 
GRB 071003 were combined into groups of three to eight 
images. The final KAIT photometry for the GRB 071003 
OA is listed in Tabled The reported error bars are the 
uncertainties in PSF-fitting photometry and those in the 
calibration process, added in quadrature. A plot of the 
KAIT photometry, along with measurements from other 
telescopes during the same timespan (with BAT data 
overplotted and fitted by a chromatic model described in 
§ 4.3) is presented in Figure [5] 

3.3. AEOS Data Reduction 

The ABC images were processed using dark subtrac- 
tion only. Because of highly variable stray light and vi- 
gnetting, we did not apply a flat field to th ese images. 
We used SExtractor (jBertin fc Arnoutslll996t ) to find all 
the sources in the images, from which we were able to 
determine the astrometry. 

We employed the NN2 flux difference method 
(jBarris et al.l 120051 : hereafter the NN2 method) for con- 
structing the AEOS light curve. The NN2 method also 
uses image subtraction to measure the fluxes for a vari- 
able source, but it does not designate one particular im- 
age as the template. Instead, given N total observations, 
the NN2 method solves for the vector of fluxes from the 
individual images using the antisymmetric matrix of flux 
differences from the N(N — l)/2 distinct possible sub- 
tractions. Compared to the template image subtraction 
method, the NN2 method takes all the available informa- 
tion from the images into account, and is less susceptible 
to possible noise associated with a single template im- 
age. To avoid a large number of image subtractions, we 
combined the original 238 AEOS observations into 39 im- 
ages. For the first 228 images, each set of six consecutive 
images is combined into one. The last 10 images are com- 
bined into a single image. We compared the results from 



the NN2 method to those from a traditional template 
image subtraction method (bottom panels of Figure [4| 
and found them to be consistent with each other. 

To calibrate the AEOS data to the standard photom- 
etry system, we used the KAIT i?-band data during the 
overlap period and assume that the unfiltered AEOS 
data have no color term to the R band. 26 The final 
AEOS photometry is listed in Table [3l The reported er- 
ror bars are only those output by the NN2 method, and 
do not include a possible large systematic error due to 
calibration. If the throughput of the AEOS telescope in 
the unfiltered mode is not drastically different from that 
of KAIT, we estimate the systematic error to be ^0.07 
mag when the GRB OA was bright (t < 20 minutes), 
and ~0.15 mag when the GRB became faint (t > 40 
minutes). The systematic errors can be much higher if 
the unfiltered throughput is very different for the two 
telescopes. 

3.4. Keck 1/ Gemini- S Data Reduction 

Due to the large aperture of the Keck I and Gemini- 
S telescopes, the bright star close to the GRB 071003 
OA produces numerous diffraction spikes, as well as two 
large blooming spikes along the readout direction. Be- 
cause the orientation, width, and intensity of the spikes 
change with the seeing conditions, the exposure duration, 
and the time of the observations, it is difficult to cleanly 
remove them using the template image subtraction or 
the NN2 method. However, due to the high resolution of 
these images, the spikes are well sampled and show dis- 
tinct axial symmetry. We developed a saturation spike 
subtraction method, in which we divide the image of the 
bright star in half, flip the right side, and subtract it 
from the left side. Due to the symmetry in the spikes, 
this subtraction process leaves a reasonably clean region 
around the GRB OA. PSF-fitting photometry was then 
performed on the GRB OA in the spike-subtracted im- 
ages, and on a series of local standard stars. The Lick 
calibration is used to calibrate the Keck I and Gemini-S 
instrumental magnitudes to the standard system. 

The final Keck I and Gemini-S photometry is reported 
in Table SI The error bars of the magnitudes are the un- 
certainties from the PSF-fitting photometry and those 
in the calibration process added in quadrature. One spe- 
cial data point is the Keck I HIRES guider image at 
t = 9523.7 s because it bridges the early KAIT/ AEOS 
data to the late-time Keck I and Gemini-S observations. 
The GRB OA was well detected in the guider image, but 
because the image has a small FOV and is unfiltered, 
photometric calibration becomes particularly difficult. 

We have used three methods to calibrate the measured 
instrumental magnitude of the OA after the guider im- 
ages were processed with the saturation spike subtrac- 
tion method: differential photometry between the AEOS 
unfiltered data and the guider images, photometric cali- 
bration to about half a dozen stars in the HIRES guider 
images using the KAIT unfiltered images, and photomet- 
ric calibration to these stars using the Keck I i?-band 
images. The measured i?-band magnitudes from these 

26 We attempted to quantify the color term of the unfiltered 
AEOS data to the standard R system using the local standard 
stars in the field of GRB 071003, but found no apparent correlation 
between the scatter of the (unfiltered — R) differences versus the 
colors of the stars. 



three methods show a scatter of ~0.25 mag, and their 
average value and uncertainty are listed in Table [H 

3.5. P60 Data Reduction 

The P60 data reduction is presented in this section be- 
cause it employs several methods (illustrated in Figure[4]) 
discussed earlier in the paper. We obtained template im- 
ages for the field after the OA of GRB 071003 has faded. 
However, the saturation spikes of the bright star close 
to the GRB ruined the template images in the R and 
i' bands, so we were only able to run image subtraction 
for the data in the g and z' bands. We also employed 
the saturation spike subtraction methods as described in 
§3.4. Although P60 does not have the resolution of the 
Keck I and Gemini-S telescopes, subtraction of half of 
the saturation spikes helped to clean up the background 
of the OA considerably. 

We also applied a third method to reduce the P60 data. 
Due to the richness of stars in the GRB 071003 field and 
the large field of view of the P60 camera (12'. 9 x 12'. 9), 
we were able to pick a star that is close in brightness 
(within 0.1 mag in all filters) and thus has similar satu- 
ration spikes to the bright star close to GRB 071003. 
The chosen star is located at a = 20' I 07 m 14 s .84, 5= 
+10°53'59".8 (equinox J2000.0), which is 136". 7 west 
and 172".0 south of the GRB 071003 OA. By slightly 
scaling the PSF of this bright star and subtracting it from 
the star close to the GRB, we were able to largely remove 
the complicated background around the GRB OA. 

PSF-fitting photometry is applied to the images after 
different ways of image subtraction, and the Lick cali- 
bration is used to calibrate the instrumental magnitudes 
into the standard system. The final photometry from 
the P60 data is listed in Table which is the average 
of the spike and bright star subtraction methods. The 
results from the template image subtraction method are 
not considered because the method can only be applied 
to a subset of filters, but they are consistent with the 
other two methods within measurement uncertainties. 

3.6. Keck AO Data Reduction 

On 2007 October 19 (starting at UT 05 :14) we observed 
the G RB 071003 OA with the NIRC2 (|Van Dam et all 
120041 ) narrow-field camera (0".01 pixel -1 ) on Keck II us- 
ing natural guide star adaptive optics (NGS AO). While 
the extremely bright nearby star greatly complicated the 
optical analysis, it was ideal to be used as the natu- 
ral guide star during NGS AO imaging. We took 15 
science exposures, each of 60 s and 2 coadds, resulting 
in a total integration time of 30 minutes. The images 
were reduced using standard techniques, including dark 
subtracting, flat fielding, and filtering for deviant pix- 
els. Each frame was dewarped using the recommended 
method for NIRC2, and the resulting images were regis- 
tered to a common origin and combined. 

The GRB OA is well detected 2 weeks after the burst, 
as shown in the final combined image in Figure [6l To 
measure the brightness of the OA, we created a model of 
the PSF using short-exposure, unsaturated images of a 
nearby Two Micron All Sky Survey (2MASS) star [K„ = 
12.011 ± 0.024 mag, d = 7.8"), taken immediately prior 
to the science exposures. We then subtracted this model 
PSF from the OA. With the same 2MASS star as the 
photometric calibrator, we measure the OA to have K' 
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Fig. 6.— NGS AO image of the GRB 071003 field taken with 
Keck II on 2007 October 19, 16 days after the burst. The FOV 
is approximately 12" X 10". The afterglow is well detected with 
K' = 21.58 ± 0.03 mag. No host-galaxy emission is detected. 

= 21.65 ± 0.10 Vega mag. (Galactic reddening of Ak< ~ 
0.05 mag is negligible along this sightline and has not 
been applied.) 

3.7. Keck LRIS Spectroscopy Reduction 

We obtained low-resolution optical spectroscopy of the 
optical afterglow of GRB 071003 on 2007 October 4.335 
using the LRIS on the Keck I telescope. A pair of 600 s 
dithered exposures was taken under clear conditions at 
airmass 1.2 with 0.6" seeing. We used both the blue 
and red arms of LRIS, with the light split by the D680 
dichroic. The 300/5000 grism on the blue side gave a 
spectral resolution of 8.4 A over the range 3300-6500 A. 
We used the 600/10000 grating to achieve 4.1 A resolu- 
tion over the range 6500-8 630 A. The spectrophotomet- 
ric standard star Feige 110 ()Stoneill977h was observed the 
following night in the same setup. Intermittent clouds 
were present the night of the standard-star observation, 
so the absolute flux scale is unreliable. 

The long, 1.0"-wide slit was oriented at a position an- 
gle of 10° for the after glow observation s, which was not 
the parallactic angle (jFilippcnkol 1 1982t h However, the 
Cass egrain Atmospheric Dispersion Compensator mod- 
ule dPhillips et alJl2006h was mounted, so the derived 
spectral shape should be reliable. The exception is in 
the spectral range of 6000-6500 A, where second-order 
blue light contamination is prominent in the spectrum 
of the standard star. An attempt was made to cor- 
rect for the contamination, but the spectral slope in 
this section is more uncertain than in the rest of the 
spectrum. We also fitted an extinction-corrected power 
law to the flux-calibrated spectrum (excluding line and 
second-order contaminated regions) in an attempt to es- 
timate the spectral slope, but the estimated slope of 
/„ oc j/ -0 - 87 differs significantly from the spectral slope 
estimated from multi-band late-time photometry (§4.7). 
This may be due to continuum contamination from the 
nearby star in the spectrum (which is difficult to properly 
remove) , so we do not further consider this spectroscopic 
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Fig. 7. — Spectrum of the GRB 071003 afterglow covering the 
full observed spectral range. The spectrum has been flux-calibrated 
and corrected for Galactic reddening of E(B — V) = 0.148 mag. 
The inset shows an expanded view of the region surrounding the 
Fe and Mg absorption system at the burst redshift. A power-law 
continuum was fitted to the regions of the spectrum shown in green, 
chosen to avoid strong absorption lines and the wavelength range 
contaminated by second-order blue light. The thick solid blue line 
shows the resultant fit (fx oc A -1 ' 13 , or /„ ex i/~ ' 87 ), but it differs 
in slope from our more reliable fit to the broadband photometry; 
thus, it is used only to normalize the spectrum. 

spectral index. 

The largely featureless spectrum (Figure [7]) has a S /N 
> 5 pixel -1 down to ~3500 A. There is no apparent 
absorption by the intergalactic medium at these wave- 
lengths, yielding an upper limit to the redshift of the 
burst of z GRB < (3500/1216) - 1 = 1.88. Numerous 
metal-line absorption lines (but no emission lines) are 
visible in the spectrum. We have fitted the equivalent 
widths of all > 5<r features in the normalized spectrum 
using a Gaussian profile and report the rest-frame values 
in Table 

We previously presented (|Perlev et al.l 12007a! ) analy- 
sis of this spectrum, identifying Mg II absorption sys- 
tems at z = 0.372 and z = 1.100. A VLT s pectrum ac- 
quired the same night (Fug azza et al.|[2~007| ) identified a 
third absorption system at z = 0.937, which is confirmed 
by our observations. These are the only strong absorp- 
tion systems in the data, and previously we considered 
it likely that the z — 1.100 system originated from the 
host galaxy (Figure [8]). Surprisingly, however, a more 
thorough investigation revealed a fourth, weak absorp- 
tion system at a higher redshift of z = 1.604 (Figure 
Contrary to our expectation, the gas at this redshift has 
the weakest Mg II absorption of the four systems. 

This is remarkable: absorption lines associated with 
GRB environments are generally very strong with rest- 
frame equivalent widths exceeding several angstroms 
(|Savaglio et all 120031 : IFrochaska et aT1l2008a|K Figured] 
also indicates, however, the presence of fine-structure 
Fe II transitions at this redshift. With the excep- 
tion of active galactic nucleus environments, these tran- 
sitions have only been identified in gas sur rounding 
the GRB phenomenon (|Prochaska et al.l 120061 ). These 
transitions are excited by the GRB afterglow itself 
throu gh indirect ultra violet pumping (jProchaska et al.l 
120061: IVreeswiik et al.l 120071 ) of gas in the interstellar 
medium (ISM) of the host galaxy. Altogether, the coin- 




CIV 1548. 1550- 



' 0.372: Fell 258d. 2600; Mgll 2796, 2803. Mgl 285. 



3500 



3600 



3700 



3800 



3900 



4000 



1.0 



0.8 - 




0.6 

o 

Z 

0.4 



:=0.937: Mgll 27%. 2803; Mgl 2852 



1.100: Fell 2586. 2600: Mgll 27%, 2803. Mgl 2852 
:0.372: Call 3934, 3969 



z GRH =1.604: 
Fell 2344, 2374, 2382 
Fell* 2396. 2405, 2411, 2414 



5400 



5600 



5800 o 6000 
Wavelength (A) 



6200 



Fig. 8. — Portions of the normalized Keck LRIS spectrum of the 
GRB 071003 afterglow. We mark the positions of several metal 
absorption-line features from four distinct extragalactic systems 
including a series of Fe II and Fe II* transitions associated with 
the host galaxy of GRB 071003 {zqrb = 1.604). Note that the 
Ca II doublet marked as Galactic may be due to the very bright 
Galactic star offset by 6.5" from GRB 071003 as opposed to the 
Galactic ISM. 

cidence of (1) the absence of any higher-redshift absorp- 
tion systems in our spectrum, (2) the positive detection 
of fine-structure Fe II transitions, and (3) the absence of 
intergalactic medium absorption at A > 3500 A estab- 
lishes z = 1.604 as the redshift of GRB 071003. 

It might seem unusual to have detected fine-structure 
Fe II transitions in such a late-time spectrum (t 
24.3 hr). Because the lines are excited by the GRB 
afterglow, they will decay as the afterglow fades on 
hour-long time s cales ( Dessauges-Zavadskv et al.l 120061 : 
IVreeswiik etaLl [20071: iD'Elia et al.l I2008D . The pres- 
ence of fine-structure transitions in our spectrum, how- 
ever, is consistent with the late-time rebrightening of 
GRB 071003 provided that the gas lies within a few kilo- 
parsecs of the GRB. In Figure[H]we present a velocity plot 
of strong resonance- line transitions for z — zgrb ■ We re- 
port the positive detections of C IV A1548, Fe II AA2382, 
2586, 2600, and Mg II A2803, and we note probable 
but statistically insignificant absorption at Al II A 1670 
and Mg II A2796. The rest-frame equivalent widths are 
among the lowest ever recorded for the ISM surround- 
ing long-duration GRBs. The equivalent width of Mg II, 
for example, is full y an order of ma gnitude below the 
general population fCenko et al.ll2008l) . with the sole ex- 
ception of GRB 070125, and the equivalent width for the 
C IV gas (Wwi8 = 0.2 2 ± Q.Q6 A), represents the lowest 
measurement to date ()Prochaska et al.ll2008af) . 

4. RESULTS AND MODELING 

4.1. Light Curve: General Observations 

The multi-color photometric evolution of the GRB 
071003 OA is shown in Figure LTOl fitted by our preferred 
model (described later). Visual inspection of the light 
curves reveals what appear to be three distinct compo- 
nents: an overall power-law decline that has already set 
in by the very first measurement at 42 s, a small "bump" 
feature at ^120-600 s, and then a dramatic, but unfor- 
tunately not well sampled, rebrightening starting around 



9 



1.0 
0.9 
0.8 

1.0 

0.9 

X 
3 

E 0.8 
-a 

_N 

13 1.0 

a 

0.8 

1.0 
0.9 
0.8 



CIV 1548 



CIV 1550 




A1II 1670 



Fell 2382 





Fell 2586 



Fell 2600 



Mgll 2796 




Mgll 2803 



-500 



500 



-500 



500 



Relative Velocity (km s ) 



Fig. 9. — Velocity plot of strong, resonance-line transitions for 
gas associated with GRB 071003 (zqrb = 1-60435). These lines 
are very weak (note the ordinate scale) with rest-frame equivalent 
widths of 100-200 mA (Table 0. Indeed, the C IV absorption 
is the weakest yet reported for a GRB afterglow (Prochaska ct al. 
I2008an . 

3000 s that dominates the remainder of the evolution. 

The bump feature appears to be additive only: fitting a 
single power law to measure the decay index (t~ a ) for the 
clear-band data both before this period and after it, the 
power-law indices (a = 1.47 and a — 1.49, respectively) 
are fully consistent with each other and with the overall 
decay index over both periods (a = 1.48). 

The rebrightening is more difficult to characterize. We 
have no observations between the Keck I HIRES guider 
point at t w 2.6 hr and our observations the second 
night; moreover, the points reported in the GCNs are 
highly discrepa nt. An opt i cal .R- band limit is reported 
at t ~ 4 hr by iShih et ail (|2007f ) , which seems to con- 
tradict the rebrightening trend suggested by the AEOS 
data and guider point. It is unlikely that the OA would 
show such a dramatic drop (> 3 mag) in a short time 
interval at such late times, so we suspect that the OA 
might be heavily contaminated by the bright nearby star 
an d was not resolve d in the Lulin 1-m telescope images 
of IShih et all (|2007| ). On the other hand, the P;band 
detection at t ~ 7.5 hr reported by Mi sra et al.l (|2007l ) 
supports a rebrightening but is several magnitudes above 
the extrapolated light curve at this time, seemingly far 
too bright to be consistent with our observations. Cali- 
bration and the contamination from the bright star are 
the likely causes of the discrepancy. 

4.2. Optical to Gamma-Ray and X-Ray Comparison 

The BAT and XRT light curves we derive for 
GRB 071003 are also shown in Figure[10] Unfortunately, 



because Swift was still in the process of returning to nor- 
mal operations after its gyro failure (Gchrcls 2007), au- 
tomatic slewing to GRB 071003 was disabled at the time 
when the GRB was detected. As a result, there were 
no prompt XRT observations for GRB 071003, leaving a 
long gap in the gamma-ray/X-ray light curve at t — 200- 
20000 s. In particular, there are no X-ray observations 
until approximately the peak of the rebrightening in the 
optical band. Nevertheless, direct comparison of the data 
available reveals three relevant facts. 

First, there is no obvious optical prompt counterpart 
to the last spike of the gamma-ray light curve. However, 
this spike is nearly contemporaneous with the much more 
slowly rising optical bump feature; we return to this pos- 
sible connection in our later modeling (§4.5). 

Second, at late times the X-ray light curve declines 
as a power law with decay index consistent with that 
observed in the optical. A simple power law fits the data 
well, with a best-fit decay index of ax — 1-68 ± 0.04. In 
addition, the late-time OA behavior (after t w 5 x 10 4 s) is 
consistent with a single power-law decay with an index of 
ao = 1.72 ±0.07, fully consistent with this value. As we 
note later, an extrapolation of the X-ray spectral index is 
also consistent with the optical observations, suggesting 
that at late times there is no need for an additional X-ray 
contribution (such as inverse Compton) or large amounts 
of host-galaxy extinction. 

Finally, while the gamma rays are scaled arbitrarily in 
Figure I10[ we note that if we extrapolate the gamma- 
ray spectrum into the X-rays to compare the BAT and 
XRT light curves, the evolution between the end of the 
prompt emission and the start of the XRT observations 
is nearly consistent with a simple extension of the late- 
time XRT power law back to earlier times, without a 
need for a rebrig htening or break. H owever, Swift has 
shown previously (|Nousek et al.l[2CK)6l ) that early-time X- 
ray light curves can conceal a wide variety of complex 
features, so we will not speculate further as to whether 
or not this was actually the case. 

4.3. Detailed Optical Modeling 

The procedure used to model the optica l light curve 
is gen erally the same as that employed by iPerlev et alJ 
(2008), but further generalized. For our fit model, we 
emp loy an unbroken power-la w decay (component 0) plus 
two iBeuermann et ail (|1999l ) functions (broken power- 
law pairs, components 1 and 2), but allow for different 
values of the functional parameters for each filter and 
component. The functional form is 



F v = F QM (t - dt )- ao 

+F 1 v (0.5 (*^L)~ s i«i,b 
+F 2 ,„ (0.5 (*^*a)- s ^2, 6 
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tpl I I 

t — dti \ — S2CX2,, 
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where for component 0, ao is the power-law decay in- 
dex, and dto is an adjustment to the Swift /BAT trigger 
time. For cmponent 1, ct\^ and ai ja are the power-law 
decay indices for the rising and declining components, re- 
spectively, dt\ is an adjustment to the Swift /BAT trigger 
time, t p i is the time of the peak flux, and si is the sharp- 
ness parameter. Component 2 has a similar function as 
component 1. 
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Fig. 10. — Multi-color, early through late-time light curves of the OA of GRB 071003. The magnitudes are offset according to their early- 
time colors, showing the color evolution between early and late times. Overplotted colored curves indicate the best-fit three-component, 
color-evolution model described in the text; the dashed lines represent the individual components that compose this model (a uniform 
power-law decay, a chromatic early-time bump, and a monochromatic late-time rebrightening). The X-ray and gamma-ray afterglows are 
also overplotted for comparison. The gamma-ray light curve is scaled arbitrarily; if scaled based on the likely gamma-to-X-ray spectral 
index it would fall on or near the extrapolation of the X-ray light curve back to early times. 



Fitting this function with no constraints generates un- 
realistic results because of non-uniform sampling in dif- 
ferent filters. However, we can make the following phys- 
ically motivated assumptions to tie specific parameters 
and produce more physically meaningful results. 

1. We assume that the temporal decay index at any 
given time is independent of the filter, as is implicit 
in the notation (a does not depend on v). This 
means that the color of a component cannot change 
except while the light curve of that component is 
breaking. 

2. Most importantly, we assume that differences be- 
tween the spectra of the various model components 
can be described by changes in the power-law index 
of the intrinsic spectrum, modified by an arbitrary, 
but fixed, extinction law. Mathematically, this con- 
straint is expressed as Fi V = v Al3ij Fj^. Physi- 



cally, this assumption requires that external effects 
such as extinction, which might cause the spectrum 
of any component to deviate from a power law, 
affect all components equally and are not time- 
dependent. The extinction law itself (as well as 
the absolute underlying index of any specific com- 
ponent) is fully general and can be fitted according 
to various models later. 

3. In addition, we assume that the rising segments of 
each component are also power laws, but not neces- 
sarily the same power laws as the falling segment, 
to allow for chromatic breaks. This imposes the 
following condition: (^) = (^)A/WAa 6a _ jjere 

b and a refer to "before" and "after" the break of a 
specific component (0, 1, or 2), where the compo- 
nent index is omitted for clarity, and x and y refer 
to two different filters. 
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Fitting is performed, under these assumptions, using 
the IDL package mpfit 27 . 

The assumptions involved in these constraints are, of 
course, oversimplifications for the full array of models 
that might be considered. In particular, this model al- 
lows only one break per component, but with an evolving 
synchrotron light curve plus a jet we may expect as many 
as three. However, it has the advantage of being simple 
and generates a single physically motivated parameter 
quantifying color change over each component. 

We perform a variety of fits under varying combina- 
tions of assumptions. Some of the possibilities we con- 
sidered include the following: 

1. Forcing the bump (component 1) to have the same 
color as the uniform decay (component 0), or al- 
lowing it to be a different color overall. 

2. Forcing the bump itself to be achromatic over its 
evolution, or allowing it to contain a chromatic 
break. 

3. Forcing the late rebrightening (component 2) to 
have the same color as the uniform decay, or al- 
lowing it to have a different color. 

4. Fixing dto for the early steep decay to be zero (the 
BAT trigger time), or allowing it to be free to vary. 

5. Fixing dt\ for the bump component to be zero, to 
be equal to the beginning of the prompt-emission 
pulse that is nearly contemporaneous with it, or 
allowing it to be free to vary. 

6. Fixing dti for the late rebrightening to be zero, or 
allowing it to be free to vary. 

The results under various combinations of these as- 
sumptions are presented in Tables [5] and [5] We discuss 
the implications of these results in the remainder of the 
paper. 

4.4. Color Change 

Detection of a GRB afterglow in filtered observations 
during prompt emission, as was the case here, is rare. 
The situation is even more intriguing since our multi- 
color prompt OA observations show an apparent bump 
feature (component 1) that is nearly contemporaneous 
with a rebrightening pulse in the gamma-ray light curve. 
Therefore, it is of great interest to attempt to measure 
the color of Component 1. By the same token, we have 
good spectral coverage of the afterglow both during the 
primary normal decay and during the fading of the dra- 
matic late rebrightening, and any color difference may 
shed light on the origin of these features. 

We tested for color differences in three places: between 
component (rapid decay) and component 1 (bump), 
between component and component 2 (rebrightening), 
and over the break of component 1 itself (since the ris- 
ing spectral index may differ from the falling spectral 
index). In all cases we find evidence for color variation, 
although in each case only at the ~ 2a level. The fading 
component of the bump is redder than the fading com- 
ponent of the uniform decay by A/3 = 0.75 ± 0.34, the 

27 \protect !ittp://cow. physics. wisc.edu/~craigm/idl/idl. html . 



bump feature is chromatic with a shift from the rising to 
falling component of A/3 = 1.11 ± 0.47, and the rebright- 
ening (for which we only have color information during 
the fading component) is also redder, by A/3 = 0.84 ± 
0.31. 

One must be somewhat cautious in interpreting these 
results — since different filters sample the data differ- 
ently, systematic errors that affect only one portion of 
the light curve can masquerade as color change. Data 
reduction for GRB 071003 was also challenging due to 
the presence of the nearby bright star, as detailed in §3. 
In addition, we note that the degree of spectral index 
shifts noted is dependent on the model. In spite of these 
considerations, however, we feel that our conclusion of 
color change is reasonably secure in each case. 

4.5. Energy Injection Times 

It is often unclear what time is most appropriate to 
use as to when fitting a power law to a GRB afterglow. 
Thanks to the extremely early-time clear-band data, it is 
possible to fit to and constrain this within a few seconds 
in the case of GRB 071003. This fit, notably, gives a 
to of exactly the trigger time (dto — —0.01 ± 3.01 s). 
The gamma-ray light curve (Figure [lj fluence is strongly 
dominated by the initial pulse, which rises sharply and 
peaks within a few seconds, so this is not necessarily 
surprising. 

So me authors jlBlake et al.l l2005t IVestrand et~aTll2005l . 
l2006t lYost et al.l 120071 ) have presented evidence of an 
optical component rising coincident with the prompt 
emission, although significantly longer lasting. We 
can analyze whether the bump component observed in 
GRB 071003 may be such a feature by determining 
whether or not it can be fitted with a pulse that rises 
abruptly, contemporaneous with the prompt emission. 
While our power-law model is somewhat simplified and 
the sampling of the rise is extremely poor, we find that 
it generally does not: the best-fit to is intermediate be- 
tween the trigger time and the time of the prompt emis- 
sion spike (^125 s) at dt\ = 60 ± 20 s. This is a model- 
independent result, although it rests mostly on one data- 
point: the initial y-band measurement, representing an 
integration from 97 to 117 s after the BAT trigger (~18 s 
before the rise of the prompt emission spike), lies 0.14 
mag above a simple power-law extrapolation from re- 
gions of the data excluding the bump, compared to a 
photometric error of only 0.03 mag. While it is possible 
to envision scenarios where a relatively slow optical rise 
might follow a gamma-ray pulse (any broadband feature 
with hard-to-soft evolution, or perhaps a late internal 
shock that later collides with and energizes the exter- 
nal shock), no model to our knowledge can explain why 
an optical flare would precede a gamma-ray pulse, so we 
take this as evidence that the two features are physically 
unconnected. 

While our sampling around the rise and peak of the 
late-time rebrightening is poor (and dominated by the 
difficult-to-calibrate AEOS and HIRES guider images), 
we can also attempt to fit the to for the rebrightening 
component. This is significantly different from t = 0, 
with a best-fit initial time of dt2 = 1245 ±311 s. (This 
is well short of its peak time of approximately 20 ks, so 
the effect on the light curve is minor.) No prompt-like 
fluctuations or other features are observed in the light 
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Fig. 11. — VLA radio light curve fitted to an unbroken 
power law. The uncertainties in the measurements have been 
increased compared to their statistical values to take into ac- 
count the effect of interstellar scintillation. Some contempo- 
raneous late-time optical points (scaled arbitrarily) are shown 
for comparison. 

curve in this region. 

4.6. Radio Modeling 

GRB 071003 is rare among Swift bursts for having a 
bright radio afterglow. We were able to successfully de- 
tect the afterglow at two frequencies and several epochs 
spanning ~2-20 days after the burst, including observa- 
tions nearly contemporaneous with our optical data. The 
data are plotted in Figure [TT1 

This GRB is not far off the Galactic plane, and the 
rad io obser v ations are a ffected by scintillation. Follow- 
ing [Walkii] (|1998l |2001[) . the afterglow is in the strong 
scattering regime for both X and C bands. An approx- 
imate modulation index (which estimates the fractional 
rms variation) is 0.4 in the C band and 0.6 in the X 
band, over a refractive timescale of ~0.5 days in the X 
band and 2 days in the C band. This is longer than any 
integration (so the error is not reduced by integration 
time) but shorter than the interval between exposures 
(so errors are uncorrelated). 

Radio data were fitted using both an unbroken power- 
law model and a singly broken power-law model. We 
attempted the fit both before including uncertainties due 
to scintillation and with an additional 40% flux error 
added to all C-band points and 60% error added to all 
X-band points. 

Without the additional flux errors, the unbroken model 
is a poor fit, with y 2 /v = 15.32/6. A single, monochro- 
matic break improves the fit dramatically (y 2 jv = 
2.45/4). This improvement is significant at 97.4% con- 
fidence. A monochromatic radio break of this nature is 
very difficult to explain physically. However, with scin- 
tillation flux errors folded into the light curve, we find 
that a simple power law is a more than adequate fit to 
the data (x 2 jv = 1.4/6), which may suggest that we have 
overestimated the degree of modulation somewhat. (This 
is to be expected: the modulation index calculated is an 
upper limit as it strictly applies only for a point source. 
The afterglow has a physical size, which quenches the 
scintillation modulation somewhat.) Therefore, as a fi- 



nal modification, we scaled down this additional error 
until the final \ 2 ' jv « 1. Properties of the temporal fits 
are given in Table [TU1 

The uncertainty due to scintillation is in any event 
too large to allow any firm conclusions about the light 
curve. However, since only refractive scintillation is ex- 
pected to be significant, the refractive timescale is much 
longer than the several-hour timescale of individual ob- 
servations, and the C-band observations were in all cases 
taken immediately after the X-band observations, we do 
consider the measurement of the radio spectral index 
(0n = —1.15 ± 0.42) to be trustworthy regardless of any 
scintillation uncertainty. 

4.7. Spectral Energy Distribution and Extragalactic 
Extinction 

If our modeling assumptions are accurate (or nearly 
so), we can use our model to calculate the SED at any 
time using a combination of all the data available, rather 
than restricting the measurement to a small subset of the 
photometry and filters, even if the data were acquired at 
very different times in the evolution of the GRB and the 
color is not constant. 

We calculate the SED at two epochs. First, we calcu- 
late the SED at t ~ 2.67 days after the burst, the time 
of our four-color Gemini-South observations. In calculat- 
ing this SED, we perform a slightly modified light-curve 
fit: we do not perform any filter transformations (e.g., to 
convert r to R), but we fix all non-SED parameters to 
that derived from the light-curve analysis. In addition, 
we add in quadrature a calibration uncertainty equal to 
5% in all filters, with a few exceptions. For z, we use 
a 15% uncertainty. For u, we use a 30% uncertainty, 
for reasons described earlier. Finally, for K' , we use a 
large extra uncertainty of 50% due to the possibility of 
a temporal break sometime between our last optical ob- 
servations and the AO observations. (However, if such a 
break is absent, then the K' observation is much more 
precise than is given on the plots.) Unfiltered observa- 
tions are not used. We also calculate an early-time SED 
during the "normal" power-law decay at t — 1000 s, us- 
ing a fit excluding late-time measurements and measure- 
ments during the (possibly chromatic) bump. Addition 
of uncertainties is as for the late-time SED. 

The resulting SEDs are plotted in Figures [12] and HH 
After removing the effects of Galactic extinction (but not 
yet considering non-Galactic extinction), both SEDs are 
a reasonable fit to a power law, providing a general con- 
firmation of our assumptions as well as indicating that 
the host or intervening galaxies do not impose a great 
deal of frequency-dependent extinction. In support of 
our analysis from the light-curve modeling, the spectral 
indices appear to differ from early to late times: /?iooos 
= 0.62 ± 0.33, while f3 2 . 6 7d = 1.25 ± 0.09. (These values 
are direct fits to the data and do not include the effects 
of the small amount of extragalactic extinction we do 
believe to be present, which we discuss shortly.) 

Unfortunately there were no early-time observations 
outside the optical band, since Swift was unable to slew 
rapidly. However, this GRB was observed nearly simul- 
taneously in X-rays, optical, and radio during the de- 
clining phase of the late rebrightening. Therefore, it 
is possible to calculate a coeval late-time spectrum at 
all wavelengths simultaneously. The values at 2.67 days 
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Fig. 12.— Optical SED of the GRB 071003 OA at 1000 s after 
the burst, fitted using the extinction constraints derived using the 
late-time SED. The intrinsic (pre-extinction) model spectrum is 
also shown. 
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Fig. 13.— Same as Figure [12] but for t = 2.67 d after the 
burst. The data (plus an X-ray normalization, not shown) have 
been fitted with an SMC-like extinction law, with the best-fit curve 
overplotted. The intrinsic (pre-extinction) model spectrum is also 
shown. 

(the same as the first optical-only SED, above, which is 
also contemporaneous with XRT observations and within 
about half a day of the first VLA observation) are given 
in Table E] and plotted in Figure [HI 

Even without considering host-galaxy extinction, the 
optical and X-ray observations are nearly consistent with 
a common spectral index: flo = 1-25 ± 0.09, 0x = 1.14 
± 0.12, and /3 x = 0.90 ± 0.03. This consistency, plus 
the fact that the optical and X-ray temporal decays are 
identical (ao = 1-72 ± 0.31, ax = 1-68 ± 0.05), argues 
that both X-ray and optical are in the same synchrotron 
regime and the spectrum across this range is a simple 
power law. We assume this throughout the remainder of 
the analysis. 

The deviations in the observed spectral index suggest 
the presence of a small amount of extragalactic extinc- 
tion. Because of the presence of numerous absorbers and 
the unusually weak nature of the highest-redshift absorp- 
tion system, however, the appropriate assumptions for 



modeling the extinction contribution are not clear. Al- 
though Mg II is not an exact tracer of the presence of 
dust, the extremely weak line absorption at the likely 
host-galaxy redshift of z = 1.604 suggests that the dust 
column at that redshift is nearly negligible. Among the 
remaining absorbers, the Mg II system at z = 0.372 is by 
far the strongest (by a factor of ^3 in equivalent width 
compared to the next strongest system at z = 1.10), 
and is likely to be the dominant contributor to any ob- 
served dust absorption. However, this is partially offset 
by the fact that dust at higher redshift is much more 
opaque (since the observed optical frequencies are in the 
rest- frame UV at z > 1), so for the moment we remain 
agnostic as to the actual redshift of the absorbing dust. 

We fit the optical spectrum simultaneously with the 
normalized X-ray flux of -FikcV = 0.036 ± 0.004 /jjy at 
2.67 days. This value has already been corrected for 
photoelectric absorption (§2.1), and X-ray absorption is 
not considered in the fit, allowing the gas-to-dust ratio 
to be independent of the amount of extinction, Ay. 

Four different extinction models were tested. In ad- 
dition to a control fit with no extinction, we fit for 
Milky Way-like, Small Magellanic Cloud (SMC)-like, 
and Large Magellanic Cloud (LMC)-like extinctio n using 
the parameterization of the iFitzpatrick &: Massal (1990) 
("FM") model, and a mo del for extinction in s tarburst 
galaxies parameterized bv ICalzetti et al.1 (j2000l ). In all 
cases the standard average value of the ratio of total-to- 
selective extinction Ry in the reference galaxy in ques- 
tion was used. (Fits with varying Ry were attempted, 
but lacking infrared or ultraviolet measurements we were 
unable to constrain this parameter.) We performed sep- 
arate fits assuming dust at z — 0.372, 1.100, and 1.604. 

Results are given in Table [T2l We find significant evi- 
dence (/-test: 96% confidence) for a small amount (Ay 
= 0.1-0.3 mag, depending on the model) of extinction 
along the light of sight. We cannot strongly constrain its 
nature; all four extinction laws, at each of the three pos- 
sible redshifts, give reasonable fits to the observations. 
The intrinsic (pre-extinction) spectral slope (3 is strongly 
constrained to be 0.94 ± 0.03, averaged across the dif- 
ferent models. This is consistent (although marginally, 
at about the 90% confidence level) with the absorption- 
corrected X-ray measurement of j3 = 1.14 ± 0.12. 

As expected, the spectrum turns over dramatically 
somewhere redward of the optical and is declining with 
decreasing frequency by the radio band. The radio re- 
sults are discussed further in §5.4. 

4.8. Photometric Limits on a Host Galaxy and 
Intervening Absorbers 

Neither our LRIS imaging nor our late-time NGS AO 
imaging show any evidence of extension or host-galaxy 
emission consistent with the afterglow position. We 
searched for emission from a host coincident with the OA 
position by smoothing and binning the PSF-subtracted 
AO image. No host emission was detected to a conser- 
vative upper limit of K' w 23 Vega mag. 

In our first-night LRIS image (when the seeing was best 
and contamination from the bright nearby star relatively 
minimized), a faint, extended source is visible slightly 
southwest of the OA. The same source is also visible in 
the AO image, clearly resolved into a faint galaxy with 
K' w 19 mag at an offset of 2.07" southwest of the OA. 
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Fig. 14. — Broadband SED at t — 2.67 days from radio 
through X-ray observations. The shaded region shows an un- 
broken extrapolation of the X-ray fit (90% confidence region), 
which is consistent with the optical measurements. The op- 
tical points are corrected for Galactic but not extragalactic 
extinction; a best-fit model for the effects of host-galaxy and 
intervening-galaxy extinction is shown (thin cyan line). The 
locations of the cooling break and peak frequency shown are 
arbitrarily chosen; the actual frequencies are not constrained 
by the available data except that both are located between 
the radio and optical bands. 

We know from the spectral analysis that there are at 
least four systems that intersect the sightline between 
the z — 1.604 GRB and Earth, including the host itself. 
Of these, the strongest candidate for association with the 
observed galaxy is clearly the z = 0.372 system, which 
both is closest and exhibits the strongest absorption sig- 
nature. (Unfortunately, we have no spectra of the galaxy 
to confirm this.) This source appears to be a small ir- 
regular galaxy, which at this redshift would be offset by 
~ 10 kpc (a reasonable distance to explain the observed 
absorption) and approximately 0.5 kpc in half-light ra- 
dius. 

No other extended sources are detected within 3" of 
the afterglow, so our upper limit rules out detection of 
both a host galaxy and any absorbing systems within this 
distance. The corresponding limit on a galaxy luminosity 
is only mild, compared to the known GRB host distribu- 
tion. At the presumptive GRB redshift of z = 1.6, any 
host galaxy is limited to a if-band absolute magnitude 
of M(K') = -22.2 Vega mag. This value falls roughly 
in the middle of the typical range of previously studied 
GRB hosts, which appear to have if-band luminosities 
on the order of 0.1 L*, and are bluer and fai nter than 
typical SCUBA galaxies (|Le Floc'h et al.ll2003l ). 

4.9. Spectroscopic Constraints on the Host Galaxy and 
Intervening Absorbers 

The very weak absorption at the host redshift in our 
spectrum suggests a lower than average H I column den- 
sity along the sightline and/or a metal-poor gas. Be- 
cause of our low spectral resolution, however, the ab- 
sorpti on is unresolved and the line profiles may be satu- 
rated ( Prochaska 2006) . We may conservatively report a 
lower limit to the column densities by assuming the weak 
limit. In this manner, we estimate Nj\tg+ > 10 126 cm~ 2 



based on the equivalent width of Mg II A2803. For a 
solar metallicity gas, this implies logiVni > 10 17 cm -2 . 
This is a conservative estimate because the gas metallic- 
ity is presumably subsolar. Nevertheless, it is unlikely 
that the gas has an H I column density matching the 
values typical of most GRBs. 

In addition to the gas associated with GRB 071003, 
the afterglow spectrum reveals three foreground Mg II 
absorbers. Two of these have moderate rest-frame equiv- 
alent widths (W2796 ~ 0.7 A), but the lowest redshift 
system exhibits a very large value (z = 0.3722, W / 2796 = 
2.5 A). The incidence of such strong Mg II absorption at 
z < 0.5 has not yet been established along quasar sight- 
lines. These absorbers are very rare at z 0.5 however, 
and the incidence is declinin g with redshift (|Nestor et al.l 
120051: iProchter et"aill2006af ). The number of absorbers 
with W2796 > 1 A per unit redshift is l{z) = 0.13 at 
z = 0.5, and the incidence of absorbers with Wxim > 2 A 
is an order of magnitude lower. This implies that one 
would need to observe of order 100 quasar sightlines to 
detect a single absorber with Wnm > 2 A at z < 0.5. 
Although these are a posteriori statistics, this analysis re- 
minds one of the apparent enha ncement of strong Mg I I 
absorbers along GRB sightlines (|Prochter et al.ll2006bl) . 
Given its low redshift, this system will be an excellent 
case to perform follow-up observations and examine the 
properties of the galaxies hosting such systems (Pollack 
et al. 2008, submitted) The bright nearby star, however, 
poses a formidable obstacle for non-AO ground-based ob- 
servations. 

4.10. Energetics 

The measured gamma-ray fluence of 5.32 (—0.67, 
+0.30) xlO" 5 erg c m" 2 (Konus, 20 keV-4 MeV: 
iGolenetskii et al.ll2007f) can be converted to an isotropic- 
equivalent total energy release in the host frame: E\ so = 
3.4 (—0.6, +0.2) xlO 53 erg — well in the upper range of 
Swift events. 

No clear jet break is observed over the course of our 
observations, in either the optical bands or the X-ray, out 
to at least 6 x 10 5 s. There is a possible monochromatic 
break in the radio bands at around 8 days (7 x 10 5 s), but 
it appears likely to be a scintillati on artifact (see §4 .6). 

Using this limit, and following iSari et "all |1999), for 
a uniform circumburst medium we can calculate the 
minimum jet opening angle and minimum collimation- 
corrected energy. Using standard values for the radiative 
efficiency (77 = 0.5) and circumburst density (n = 3.0 
cm~ 3 ) (the end result is nearly insensitive to these pa- 
rameters), we have 

0» = 6.5°(^) 3 / 8 (-^) 1 /8(i±i)-3/8 ( 

d 3 cm d 2 \\j a6 erg 

(2) 

However, as we discuss later, the late-time afterglow be- 
havior in this case fav ors a wind model. Thus, following 
iLi fc Chevalier] (|2003l ) we have 

9 ict = 5.4°(^) 1 /4 ( ^ ) i/4 ( i±f } -i/4 ( Jg/H-)-V* 
d 2 erg 

The upper limit on ij et of 7 days gives a limit on the 
opening of at least 3.1 (A^/O.l) 1 / 4 deg. (As discussed 
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later in §5.4, we estimate A* « 0.1 from the broadband 
spectrum.) Therefore the collimation-corrected energy is 
at least E 7 > 2 xlO 50 (A^/0.1) 1 / 2 erg. 

It is also possible that the jet break is hidden by the 
complicated evolution of the burst, including the re- 
brightening, which would imply more modest energetics 
for this burst. However, as the late-time slope is still 
relatively shallow (a = 1.72; generally we expect a > 2 
after a jet break) we consider this relatively unlikely. 

5. DISCUSSION 
5.1. Initial Power-Law Decline 

We first turn our attention to the rapidly declining 
power law. The temporal behavior of this feature is quite 
simple, with a decay constant a = 1.466 ± 0.006 and 
no evident substructure before the "bump" or after it. 
There is no evidence of a rising component or any early 
break. The observed spectral index (3 = 0.62 ± 0.33, 
although if the extinction measured at late times is also 
present at early times (as we expect), the intrinsic index 
is actually shallower; correcting this using our preferred 
extinction model, we derive f3 = 0.29 ± 0.49. 

Especially when the decay is observed to flatten later, 
very early-time decay of this nature is often interpreted 
as a reverse shock. This seems possible — the spectral 
and temporal indices are within the range of predictions 
fo r reverse- shock m odels (specifically, the thick-shell case 
of iKobava shi 2000). However, a forward-shock origin is 
also consistent. Examining the standard c losure rela- 
tions between a and (3 (as in, for example, iPrice et all 
2002), all environment models (ISM, wind, and jet) are 
consistent with the constraints derived from the data, 
largely because the early-time constraint on (3 is poor. 
(We discuss the forward vs. reverse-shock models for this 
emission again in §5.3, in connection with the late-time 
rebrightening.) 

5.2. The Bump: Internal Shock Origin Without a 
Prompt- Emission Connection 

The bump feature is of considerable interest, since it is 
nearly simultaneous with a prompt-emission pulse. How- 
ever, as discussed earlier, the temporal analysis seems to 
disfavor the interpretation as a prompt reverberation: 
the bump seems to be already rising even before the 
prompt spike. 

Another possible explanation for the origin of this 
feature is a large density variation in the surrounding 
medium (a large clump or other discrete physical fea- 
ture in the path of the expanding shock) . The observed 
pulse width At/t m 1 is consistent with a density varia- 
tion, and the general appearance of the light curve over 
this region is reminiscent of simulations of a GRB for- 
ward shock intersecting ISM dens ity enhancements (e.g., 
Figure 3 of iNakar &; Piranl [2003). However, our obser- 
vation of possible color change across the bump would 
(if real) disfavor this hypothesis, at least in the sim- 
plest models: density variation will not change the intrin- 
sic spectrum, unless either the microphysical parameters 
or cooling frequency suddenly and significantly change. 
We consider this unlikely, although som e authors (e.g., 
lYost et al.|[200l iGranot & Kumarll2006fl have discussed 
the role of variable microphysics in previous GRB after- 
glows. 



Alternatively, the observation that the fast-declining 
component seems completely unaffected by the after- 
glow (the temporal indices before and after are effectively 
identical) leads us to interpret the bump as originating 
from a distinct emission episode — given the rapid rise 
and fall and the hint of blue-to-red evolution we sug- 
gest that it arises from internal-shock emission. Hard- 
to-soft evolution and an underlying power-law decay not 
affected by the flare have also been seen in X-ray flares 
(|Butler fc Kocevskill2007at IChincarini et al.ll2007f ). We 
also note that earlier studies of GRB prompt emission 
have shown pulses observed at lower energy to be br oader 
than those at higher energy (jFenimore et al.iri995l ): this 
trend may continue into the optical band. The broader, 
smoother profile of this pulse relative to the much faster- 
evolving X-ray flares may in this case illustrate important 
attributes of the emission — either from viewing effects 
or resulting from the physics of the emission itself. 

5.3. The Late Rebrightening 

The rebrightening phase of this burst is quite dramatic. 
While our observations do not sample the peak of the 
emission, a fit with a reasonable assumption of the sharp- 
ness parameter suggests that the flux increased by ap- 
proximately 1 mag, and the amount of integrated optical 
flux released during the rebrightening is comparable to 
or more than that emitted by the early afterglow. A rise 
in optical flux of more than a magnitude at intermediate 
times (well after the end of prompt emission, but before 
any supernova component) has to our knowledge been 
seen in only a handful of previous ca ses: GRBs 97050 8 

, 041219 A dBlake et al.l2005[) . 
, 07042 (jJelmek et al.ll2007fi . 
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and 070311 (|Guidorzi et al.ll200l . 

The rebrightening is also notable because it appears 
to differ subtly from the early decay, even though the 
evolution of both curves is generally quite simple. The 
decay index and spectral index both steepen, by Aa = 
0.25 ± 0.14 and A/3 = 0.80 ± 0.30, respectively. Assum- 
ing a synchrotron spectrum, there are only two possible 
origins for this — the optical band is in different syn- 
chrotron regimes at different times (specifically, v < v c 
before cooling, and v > v c after cooling, consistent with 
the changes observed), or because of a shift in the elec- 
tron index p by approximately Ap = +0.4. 

We consider several physical origins for the rebrighten- 
ing feature: the appearance of the forward shock when 
the burst ejecta first decelerate against the ISM, the 
late-time peak of a pre-existing forward shock due to 
evolution of the critical frequencies, impact of the for- 
ward shock through a density variation, and rebrighten- 
ing caused by a refreshed shock. 

Appearance of forward shock — When the GRB ejecta 
first begin to sweep up an amount of matter from the 
ISM comparable to the energy in the ejecta, they begin 
to decelerate, and reverse and forward shocks are propa- 
gated back into the ejecta and forward into the ISM, re- 
spectively; depending on the Lorentz factor, both shocks 
can then rise very quickly. We consider this scenario 
extremely unlikely to be relevant, since by necessity the 
forward and reverse shocks must rise simultaneously, and 
there is no explanation for the bright early-time compo- 
nent in the burst — save for a prompt model connected 
with internal shocks, but as we have already shown, there 
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is no evidence linking the early optical behavior with the 
high-energy emission. 

Spectral peak of existing forward shock — A more rea- 
sonable model postulates that the reverse and forward 
shocks both formed extremely early, but because they 
evolve differently (the reverse shock, whose synchrotron 
parameters are boosted down by factors of 7 2 , begins 
to fade immediately while the forward shock will rise 
at lower frequencies), the reverse shock fades rapidly, 
while the forward shock can rise and peak when the syn- 
chrotron frequency v m passes through the optical band. 
This model has, for example, been invoked to explain 
early-time bumps in the light curves of GRB 021004 (e.g., 
Kobavashi fc Zhand [2003) . GRB 050525A dShao fc Pail 
2005D, and GRB 080319B (jBloom et all l2008f) . which 
level off significantly (but do not rebrighten) at around 
10 4 s. However, this model is problematic here: although 
we have only sparse observations of the rebrightening, 
the observed rising temporal index of a = —1.12 ± 0.16 
is far too fast to be consistent with a rising phase of a 
forward adiabatic shock, which predicts F cx i( 2_s )A 4 ~ s ) 
(= t 1 !" 1 for a constant-density ISM and constant for a 
wind). Therefore, the synchrotron peak of the forward 
shock alone cannot explain this feature. 

Density variation — A third possibility, not invoking 
the transition between reverse and forward shocks, might 
be a dramatic density variation: for example, the impact 
of the shock wave into a previously ejected circumstellar 
shell, or emergence of the shock from a low-density cav- 
ity into a dense external medium. Density fluctuations 
have been successfully invoked to ex plain low-level varia - 
tions in several previous studies (e.g. JLazzati et aLll2002D 
and the timescale of the rebrightening (A t/t ~ 1) is con- 
sisten t with a density- fluctuation origin (Nak ar fc Piranl 
I2003T) - However, in this case we would expect neither a 
change in the spectral index (as is probably observed) 
nor such a slow decline after the peak, with a tempo- 
ral index that differs significantly but only slightly from 
the value of the initial decay. Furthe rmore, detailed nu- 
merical studies by numerous authors (iHuang et al1l2006l : 
INakar fc Pirar][200l iNakar fc Granotl 120071 ) have failed 
to reproduce anything but the smallest rebrightening sig- 
natures in previous GRBs using density variations. 

Multi- component jet — The com plicated light curve of 
GRB 030329 has been interpreted (jBerger et al.ll2003[ ) as 
the result of two separate forward shocks, arising from 
two different jet components: a narrow, highly rclativis- 
tic jet whose emission peaks extremely early, plus a wide, 
more mildly relativistic jet that dominates the late-time 
and radio evolution. Could this model conceivably ex- 
plain the observations of GRB 071003? While a complete 
analysis is beyond the scope of this paper, we note that 
the observations do seem consistent: the similarities of 
late-time decay of both rapid and late-time components 
are naturally explained, the timescale of our rebright- 
ening is similar to that observed in GRB 030329, and 
(notably) the most significant criticism of the two-jet in- 
terpretation of GRB 030329 (that th e rebrightening rose 
too rapidly and peaked too sharply — iHuang et al.ll2006l) 
does not apply here: the rebrightening in this case is 
much smoother than that observed for GRB 030329. 

Refreshed shock — Finally, we consider the possibility 
that this feature is due to a discrete energy reinjection 



energizing the forward shock, such as via a slow-moving 
shell that catches up to the forward shock at late times 
after it decelerates. This seems consistent with all obser- 
vations, although largely by virtue of not making strong 
predictions; by invoking a customized pattern of energy 
reinjection at the right times, a very broad space o f light 
curve behavior can be modeled (jHuang et al.ll2006T) . We 
do note that a large, sudden rebrightening of this na- 
ture may also produce a (second) reverse shock, which 
would be observable in radio and decline rapidly with 
time. The radio flux does in fact decay somewhat (in 
contrast to the expectation from a forward-shock model, 
where the radio flux is constant or rising), and the mea- 
sured a — 0.33 ± 0.10 is not far from the predicted decay 
constant for a rev erse shock of a » 1/2 in the v < v m 
frequency regime (Kobayashi 2000). However, the radio 
decay could conceivably be due to other effects (e.g., late 
jet break), and without an independent measurement of 
the synchrotron peak frequency v m and late-time Lorentz 
factor r we are unable to further constrain the presence 
or absence of such a feature with the limited observations 
available. 

We therefore find that only the multi-component jet 
and refreshed shock models are consistent with all avail- 
able data. Unfortunately, we do not have sufficient ob- 
servations during the rising phase of the rebrightening 
to distinguish the two models; in particular, we can set 
no constraints on the color evolution and lack a detailed 
light curve of the rise to peak of the rebrightening. We 
do note that the X-ray observations are already decaying 
well before the (probable) optical peak by an extrapola- 
tion of our observations (Figure I10[) , which may suggest 
hard-to-soft evolution in this feature as well. However, 
as noted earlier, the X-ray decay extrapolates back to 
the BAT light curve without explicit need for a rebright- 
ening, so without earlier X-ray measurements this asso- 
ciation is speculative. 

5.4. Environmental Constraints 

In the simplest models, the late-time light curve of any 
GRB is fixed by a number of basic parameters: micro- 
physical parameters eb (the fraction of energy in mag- 
netic fields), e e (the fraction of energy in electrons), and 
p (the electron energy index); macroscopic parameters 
Ek (the blastwave energy) and 9j (the jet opening an- 
gle) ; and a parameter quantifying the density of the sur- 
rounding medium, n (for a uniform density) or A* (for 
an r~ 2 density profile). Our broadband observations 
(spanning from radio to X-rays) should, in principle, al- 
low us to firmly constrain most of these parameters for 
GRB 071003 — or, more accurately, to its late rebright- 
ening phase, as this component is dominant at late times. 

The indices a and (3 are both well constrained at late 
times in the optical through X-ray bands, thanks to the 
wide range of temporal and spatial sampling: exo+x = 
1.71 ± 0.14, Pox = 0.93 ± 0.03. Two environment mod- 
els satisfy these constraints within 90% confidence: a 
wind-driven medium (p oc r~ 2 ) in which p « 2.9, and a 
model in which the jet break has already occurred with 
p fa 1.9 (but consistent with p = 2). Notably, ISM mod- 
els are a poor fit: the late-time decay rate is too fast 
for the shallow spectral index. The radio observations 
appear to support this conclusion: the rising light curve 
predicted by the ISM model is clearly ruled out, and 
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while the slow radio decay (<xr = 0.33 ± 0.1) is inconsis- 
tent in detail with the wind prediction of constant evo- 
lution as well, it is conceivable that variations from an 
exact s = — 2 profile, an additional source of radio emis- 
sion at early times (e.g., a reverse shock), or a soft jet 
break at t rj 5 days may explain this difficulty. 

The apparent spectral index of (3 w —1.1 observed in 
the radio is notable. A synchrotron spectrum is expected 
to have a self-absorbed (3 = —2 spectrum below the self- 
absorption frequency v a and a spectrum of (3 = —0.5 
above it. The fact that the observed spectral index is 
intermediate between these values and consistent with 
neither (to ^90% confidence) tells us that, if the spec- 
trum is really synchrotron, the absorption break is likely 
to be very close to these frequencies, although exact con- 
straints are difficult with only two frequencies since the 
break is likely to be quite soft. The radio evolution ap- 
pears nearly achromatic, which would argue against this 
interpretation, but considering the relatively narrow time 
and frequency window of the observations and unknown 
break sharpness, we feel that this is not a major concern. 

Because the ISM model is notably discrepant with the 
measured values of a and /3, we unfortunately cannot use 
the afterglow as a probe of the ambient density. If the 
wind model, which is more consistent with the observa- 
tions in this case, is correct, we can calcu late the param- 
eter A * using (for example) equation 2 in lChevalier fc Ll 
(EHi): 



F Vm = 20 mJy(- 



'5403 Mpc A ; V 0.1 ; 52 d 

While we have no direct measurement of F Vm , it is con- 
strained by the radio and optical observations (see Figure 
IT4|) to be ~1 mJy (within a factor of ~3). We therefore 
measure A* — 0.07(es/0.1) 1 / 2 , an interestingly low value 
regardless of the value of cb- While es is not strongly 
constrained, the absence of a cooling break between the 
X-ray and optical bands during the first 5 days (the 
cooling frequency v c increases in a wind model) requires 
e B > 0.3. 

It is possible that the wind model is inappropriate and 
the rapid optical decay is due to a jet that broke be- 
fore our multicolor late-time observations. (One possible 
criticism of the wind model is that in this case, the color 
transition between early and late times is hard to ex- 
plain; because the cooling frequency rises with time, if 
v > v c late it must have been early as well under stan- 
dard synchrotron evolution. However, because the re- 
brightening appears to be either a separate phenomenon 
or a large energy impulse that could conceivably have 
"reset" the synchrotron parameters [including v c ] to new 
values, this may not be a major concern.) No jet break 
is observed in the light curve, but it is possible that a jet 
signature was concealed by the rebrightening. This case 
would certainly rule out the wide-angle jet interpretation 
of the secondary peak and would significantly reduce the 
energetics. 

5.5. Spectral Implications on the Environment and Host 

Galaxy 

The late-time spectroscopy and imaging tell a 
coherent story: unlike the vast majority of GRBs 



(|Wainwright et all 120071 : iProchaska et all l2008al) . 
GRB 071003 did not occur in a gas-rich 28 galaxy. The 
environment is more consistent with a progenitor located 
in an outer galactic halo, or in an extremely small (even 
compared to "normal" long-duration GRB hosts) and 
gas-poor galaxy. While the possibility of line saturation 
prevents us from setting definitive upper limits, the 
column density through any host is consistent with 
being 3 orders of magnitude below typical GRB-derived 
values, and the contrast to the overall GRB population 
- which is dominated by subluminous galaxie s to b egin 
with (e.g.. iFruchter et all [20CB iFvnbo et all 12008( 1. is 
dramatic. 

While it is well established that long-duration GRBs 
generally originate from massive stars, we should be care- 
ful to ensure that our prior experience docs not blind us 
to the existence of rarer subclasses of events. We note 
that one other GRB on record, GRB 070125, had very 
similar properties: extremely low Mg II absorption and 
no coincident host (|Cenko et alj 120081 1 . as well as a very 
brigh t afterglow and extr eme energetics (E 1 = 3 x 10 52 
erg; I Chandra et al.ll2008ft. and e ven a (mild) late-time 
rebrightening ( Updike et al.ll200c? l . Both are also among 
the few Swift bursts detected at radio wavelengths. 

However, GRB 070125 and GRB 071003 show evidence 
from their broadband light curves of origins typical of 
ordinary long GRBs. In the case of GRB 070125, a con- 
stant but very high circumstellar density suggested that 
it occurred in what was locally a dense environment, not 
an empty galactic halo, despite the near absence of a 
large-scale gas signature in the spectrum. In our case, 
for GRB 071003, we find evidence of a wind-like strati- 
fied environment, a characteristic of a massive star. To- 
gether, these events appear to suggest an origin for these 
"halo" bursts similar to those of all other GRBs. 

If GRB 071003 did occur in a star-forming region, 
then there are two possibilities consistent with the ex- 
tremely small metal absorption in the spectrum. First, 
the burst may simply have formed in an extremely sublu- 
minous galaxy — necessarily, the number or distribution 
of such objects at very high redshift is not observation- 
ally constrained, but most simulations predict an abun- 
dance of small, highly sub-Galactic halos in the universe 
that could very well harbor limited star formation. Al- 
ternatively, GRB 071003 may have occurred in a tidally 
stripped tail from another, larger galaxy. In this case, 
further follow-up observations should reveal a disturbed, 
star- forming host in the close vicinity of the burst. 

Either scenario seems plausible to explain the con- 
straints derived on the burst environment. In cither case, 
if GRBs are shown to be reasonable tracers of star for- 
mation at high redshift, then future large-sample GRB 
spectroscopy missions may be able to place important 
constraints on the star-formation history of the universe 
not possible by any other means. While the sample 
size of such low-column-density GRBs is now small (two 

28 Since our measurement is based on magnesium, we are directly 
measuring the metal column, not the gas column. An alternate 
possibility, therefore, is that the host is "normal" but extraordi- 
narily metal-poor, less than 10 -2 of the average solar abudance. 
However, we consider a highly subluminous host a more likely pos- 
sibility. Both effects may be in play: low-luminosity galaxies, and 
those with low equivalen t widths, tend to be relatively metal-poor 
UProchaska et al.ll2008bh . 
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events, with GRB 061021 [Thoene et alll2006| possibly 
constituting a third example), these results are already 
suggestive that this fraction may be significant (on the 
order of a few percent), and systematic rapid afterglow 
spectroscopy should continue to increase the number con- 
siderably over the years and decades to come. It would 
be an interesting discovery if the distribution of Mg II 
equivalent widths turns out to be bimodal. 

On a related note, the existence of GRB 071003 and 
GRB 070125 may have important implications regarding 
the escape fraction of ionizing photons and the rcioniza- 
tion history of the universe. Although the relatively low 
rcdshift of these systems keeps the Lyman-a and Lyman- 
break absorption features out of our spectral range and 
preven ts us from measur ing the H I column density di- 
rectly (|Chen et al.ll2007l ). these GRBs provide evidence 
that massive stars can form well outside of gas-dense 
hosts, where there is little to shield the intergalactic 
medium from their ionizing UV radiation. If the frac- 
tion of these events is more than a few percent at z > 7, 
then such "halo" stars may in fact be primarily responsi- 
ble for the reionization of the universe. Observationally, 
spectroscopy of such events at these high redshifts may 
allow ac curate measurement o f the neutral gas fraction 
in (e.g.. iMcQuinn et al.ll2007l ) without the interference 
of saturated line profiles originating from the host galaxy. 

6. CONCLUSIONS 

Although the temporal evolution of the optical after- 
glow of GRB 071003 is complicated, our early through 
late-time photometric follow-up data clearly resolve the 
optical light curve into separate components. Observa- 
tions from KAIT during the prompt phase of the GRB 
revealed a slowly rising, slowly falling bump or flare com- 
ponent, superimposed on a simple fading power law that 
has no observable correlation with the prompt emission, 
suggesting that while early internal-shock flares can be 
observed in the optical, they are not necessarily the same 
as those producing the high-energy signatures. Our late- 
time observations revealed one of the most dramatic late 
rebrightenings ever recorded in a GRB light curve, and 
suggest that this feature is not due to a reverse-forward 
shock transition or density variation, requiring either an- 
gular jet structure or very discrete late-time re-energizing 
of the optical afterglow. This may have important im- 
plications for the interpretation of other, less dramatic 
bumps and rebrightenings at similar timescales that ap- 
pear to be common features in GRB afterglows. 

The spectroscopic study of GRB 071003 offers a cau- 
tionary tale about the standard use of Mg II to infer a 
redshift: while it is common practice to use the highest- 
redshift Mg II system observed (especially in the cases 
when the absorption is quite strong) under the assump- 
tion that the GRB host system should show significant 
metal absorption, here we have a clear case where this 
assumption is fundamentally flawed. Were the S/N of 
the spectrum worse, or the host-galaxy absorption even 



weaker by a factor of only 2-3, it is likely that we would 
have missed the higher-redshift system entirely and pro- 
ceeded with the assumption that this burst was at a red- 
shift of 1.100 instead of 1.604. In light of this fact, previ- 
ous and future GRB redshift claims based solely on iden- 
tification of Mg II absorption should be regarded with 
increased skepticism. 

The intervening absorption systems are nevertheless 
also remarkable. With three completely independent 
Mg II systems along the line of sight, GRB 071003 is 
among the most dramatic examples yet of the bizarre 
overabundance of these systems in GRB afterglows rela- 
tive to those of quasars. Further study of this sightlinc, 
especially using AO systems, may help shed light on this 
mysterious result. 
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TABLE 1 

Photometry of comparison stars in the field of GRB 071003 



ID 


O!J2000" 


*J2000" 


B 




Mb 


V 




IV 


H 






1 




N / 


1 


oUl 


yuoo 


10.9743 


16.7601 


014' 


1 4 


15.7491 


012 


1 4 


15.235( 


008' 


i 4 


14.7201 


007 


1 4 


2 


oUl 


onno 
.yuuy 


10.9357 


13.7601 


'009* 


I 9 


12.7891 


007' 


1 8 


12.251( 


'008* 


i 9 


11.7151 


009' 


1 8 


3 


OUl 




10.9249 


14.8651 


'008* 


l 9 


13.2431 


006' 


1 8 


12.3371 


008* 


i 9 


11.4511 


008' 


1 7 


4 


oUl 


.oyuo 


10.9968 


15.6861 


! 006 ! 


I 6 


14.7481 


009' 


1 8 


14.255( 


009* 


i 9 


13.8181 


008' 


) 8 


5 


OUl 


ooll 


10.9343 


15.0261 


'009' 


I 8 


13.5141 


'008' 


1 9 


12.6791 


'007* 


i 10 


11.8821 


006' 


1 7 


6 


oUl 


070Q 

.o ( yy 


10.9803 


17.0881 


009' 


) 7 


16.0501 


009' 


1 8 


15.4991 


008* 


i 8 


14.9391 


oio' 


1 6 


7 


301 


8727 


10.9562 


17.1841 


'014 ( 


) 4 


16.3721 


009' 


) 8 


15.898( 


'008* 


i 8 


15.4791 


oio' 


1 7 


8 


301 


.8705 


11.0005 


16.1241 


'oio' 


I 5 


14.7871 


'012' 


1 6 


14.044( 


'008* 


i 7 


13.3531 


'008' 


1 6 


9 


301 


.8693 


10.9253 


16.4731 


'009' 


I 8 


15.63K 


'008' 


1 10 


15.149( 


006 


i 10 


14.714( 


'008' 


1 9 


10 


301 


.8688 


10.9989 


15.1391 


'oio' 


l 5 


14.5231 


009' 


1 6 


14.1531 


'008 


i 7 


13.8041 


'008' 


1 6 


11 


301 


.8671 


10.9417 


15.3101 


009' 


I 8 


14.5371 


008' 


1 10 


14.0951 


007' 


i 10 


13.6921 


008' 


1 10 


12 


301 


.8570 


11.0033 


15.8281 


'012' 


I 5 


15.0041 


'012' 


) 4 


14.5691 


009 


i 6 


14.1581 


'010' 


1 4 


13 


301 


.8467 


10.9749 


13.4701 


'008' 


l 7 


12.8991 


'008' 


1 9 


12.554( 


'008 


i 10 


12.2201 


'007' 


1 8 


14 


301 


.8396 


10.9967 


16.4681 


'009' 


) 7 


15.7521 


'008' 


) 9 


15.33K 


'007 


i 10 


14.9361 


'009' 


1 9 


15 


301 


.8386 


10.9736 


16.7781 


009' 


I 8 


15.89K 


009' 


1 6 


15.3671 


008' 


i 9 


14.8421 


008' 


1 7 


16 


301 


.8383 


10.9800 


15.3051 


'oio' 


l 7 


14.0341 


'008' 


1 9 


13.3571 


'006 


i 10 


12.7731 


'008' 


1 9 


17 


301 


.8354 


11.0043 


17.4551 


'017' 


) 3 


16.5601 


'004' 


) 3 


15.9701 


'010 


I 5 


15.4691 


'008' 


1 5 


18 


301 


.8342 


10.9820 


16.1091 


008' 


l 9 


14.9921 


009' 


1 7 


14.3851 


007' 


i 10 


13.8321 


010' 


1 8 


19 


301 


.8334 


10.9963 


16.4721 


'012' 


) 6 


15.4301 


'010' 


1 7 


14.7861 


'008' 


i 10 


14.2541 


'008' 


1 9 


20 


301 


.8311 


11.0048 


16.227(012' 


) 4 


15.529(009' 


) 3 


15.056(007' 


I 6 


14.631(006' 


1 5 


21 


301.8279 


10.9884 


17.045 


(008 


) 7 


16.107 


(008 


) 9 


15.546 


(008 


) 10 


15.007 


(009 


) 8 


22 


301.8247 


10.9781 


14.276 


(008 


) 9 


13.825 


(008 


) 9 


13.510 


(008 


) 10 


13.223(008 


) 9 


23 


301.8221 


10.9256 


16.606 


(009 


) 9 


16.002 


(012 


) 6 


15.623 


(007 


) 9 


15.204(010 


) 7 



Note. — Uncertainties (standard deviation of the mean) are indicated in parentheses. This table has 
been truncated; additional standard stars are available in the online material. 
a In degrees. 



TABLE 2 
KAIT photometry of GRB 071003 



tstart" 


h^xp. time 


Mag 


terror 


filter 


(s) 


(s) 








42.0 


5 


12.791 


0.019 


R b 


49.0 


5 


12.999 


0.024 


R 


55.0 


5 


13.193 


0.021 


R 


61.0 


5 


13.321 


0.024 


R 


67.0 


5 


13.500 


0.019 


R 


97.0 


20 


14.465 


0.027 


V 


128.0 


20 


13.919 


0.032 


1 


157.0 


20 


14.382 


0.031 


R 


188.0 


20 


14.916 


0.034 


V 


219.0 


20 


14.121 


0.022 


1 


249.0 


20 


14.750 


0.023 


R 


279.0 


20 


15.409 


0.035 


V 


310.0 


20 


14.578 


0.030 


1 


340.0 


20 


15.401 


0.024 


R 


370.0 


20 


16.034 


0.107 


V 


401.0 


20 


15.478 


0.063 


1 


431.0 


20 


16.239 


0.082 


R 


462.0 


20 


16.853 


0.200 


V 


492.0 


20 


15.977 


0.069 


1 


522.0 


20 


16.749 


0.091 


R 


565.0 


20 


16.255 


0.077 


I 


595.0 


20 


16.849 


0.097 


R 


624.0 


20 


16.364 


0.106 


I 


654.0 


20 


17.041 


0.089 


R 


749.09 c 


3x20 


16.830 


0.113 


I 


787.68 c 


3x20 


17.362 


0.121 


R 


1007.74 c 


6x20 


17.314 


0.148 


I 


1009. 95 c 


5x20 


17.711 


0.147 


R 


1422.34 c 


8x20 


18.103 


0.154 


R 


1464.58 c 


7x20 


17.473 


0.135 


I 



a The start time of the exposure, in seconds after 
the BAT trigger. 

b The R-band photometry is derived from unfil- 
tered observations. 

c The time (s) at the middle point of several 
combined images. 
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TABLE 3 
AEOS i?-BAND 
PHOTOMETRY FROM 
UNFILTERED 
OBSERVATIONS a 



<mid b 


R 




(s) 










568.6 


16. 


.708 





.016 


681 3 


17 


.046 


0. 


.018 


794^0 


17. 


337 


0. 


020 


906 8 


17. 


.573 





.020 


1019.6 


17. 


766 


0. 


020 


1132.3 


17. 


940 


0. 


.020 


1245.1 


18. 


101 


0. 


020 


1357 9 


18. 


229 


0. 


.023 


1470.7 


18 


.339 





.025 


1583 5 


18 


.454 


0. 


.026 


1696.3 


18. 


.545 


0. 


028 


1809.1 


18. 


640 


0. 


.034 


1922.0 


18. 


724 


0. 


.033 


2034.8 


18. 


814 


0. 


.040 


2147 6 


18 


865 


0. 


037 


2260.3 


18 


.924 





.042 


2373.1 


18 


.941 





.048 


2485.9 


18. 


.941 





.064 


2598 7 


19 


.042 





.044 


2711 4 


19. 


087 


0. 


049 


2824.2 


19. 


104 


0. 


.050 


2937.0 


19. 


109 





054 


3049.7 


19. 


142 


0. 


.054 


3162.5 


19 


150 





052 


3275.3 


19 


160 


0. 


.056 


3388.2 


19 


169 


0. 


055 


3501.0 


19 


.158 





.050 


3613.8 


19 


.139 


0. 


.056 


3726.5 


19 


.185 


0. 


.055 


3839.3 


19 


206 


0. 


.056 


4062.4 


19 


.188 


0. 


.056 


4175.3 


19 


.196 


0. 


.055 


4288.1 


19 


194 


0. 


.057 


4401.0 


19. 


208 


0. 


.051 


4513.8 


19. 


.167 


0. 


054 


4626.6 


19 


.150 





.053 


4739.4 


19 


.163 





.065 


4852.2 


19 


160 


0. 


.057 


5002.6 


19. 


145 


0. 


050 



a The original data set 
was grouped and com- 
bined into 39 images. 
b The time at the middle 
point of several combined 
images. 



TABLE 4 

Keck/Gemini-S PHOTOMETRY OF GRB 071003 



^mid 

(s) 


Jiixp. time 

(s) 


Mag 


terror 


Filter 


telescope 


9523.7 


2x20 


18.59 b 


0.25 


Cr 


Keck I 


76891.8 


300 


20.32 


0.07 


Keck I 


77044.0 


300 


19.43 


0.06 


R 


Keck I 


231174.7 


450 


22.33 


0.20 


9 


Gemini-S 


231802.3 


450 


21.57 


0.32 


i 


Gemini-S 


232430.8 


450 


21.97 


0.22 


r 


Gemini-S 


233056.5 


450 


21.35 


0.40 


z 


Gcmini-S 


515855.0 


1485 


22.61 


0.30 


R 


Keck I 


516250.5 


975 


23.56 


0.30 


a 


Keck I 


517510.4 


720 


23.56 


0.45 


u 


Keck I 


604978.0 


660 


23.06 


0.50 


R 


Keck I 


605144.4 


780 


24.05 


0.40 




Keck I 


682211.1 


330 


23.56 


0.40 




Keck I 


682211.9 


660 


24.42 


0.50 


u 


Keck I 


682940.0 


840 


23.40 


0.60 


R 


Keck I 


1373589. 


1800 


21.58 


0.03 


K' 


Keck II 



a The time (s) at the middle point of the observations. 

b Measured from unfiltered images from the Keck I HIRES 

guider. 
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TABLE 5 
P60 photometry of GRB 071003 



— 7 * — 

Estart 


i^xp. time 


Mag 


terror 


Filter 


(s) 


(s) 






176.0 


60 


14.57 


0.06 


R 


261.0 


60 


15.08 


0.05 


i> 


347.0 


60 


15.37 


0.07 


z' 


432.0 


60 


16.03 


0.07 


R 


518.0 


60 


16.41 


0.07 


i' 


603.0 


60 


16.43 


0.08 


z' 


689.0 


60 


16.88 


0.09 


R 


775.0 


60 


17.12 


0.10 


i' 


860.0 


60 


17.02 


0.11 


z' 


1309.0 


120 


17.92 


0.20 


z' 


1454.0 


120 


18.71 


0.10 




1891.0 


120 


18.67 


0.28 


I' 


2037.0 


120 


19.33 


0.15 


9 


2618.0 


120 


19.44 


0.20 


9 



a The start time of the exposure, in seconds 
after the BAT trigger. 



TABLE 6 
Radio observations of GRB 071003 



0T"DaEc5 

Observation 

-2W7 

2007 
2007 
2007 

2007 Oct. 14 
2007 
2007 
2007 
2007 



frequency 
GHz 



flux density 



frror 



^mid 

hr 



Oct. 05, 

Oct. 07, 

Oct. 07, 

Oct. 12, 



3.38 
3.92 
1.03 
14.84 
23.58 
0.04 
0.01 
0.18 



Oct. 24, 
Oct. 25, 
Nov. 05, 
Nov. 07, 



42.168 


8.46 


393 


55 


91.698 


8.46 


430 


50 


92.238 


4.86 


220 


54 


209.248 


8.46 


431 


51 


271.158 


8.46 


332 


67 


519.898 


8.46 


260 


42 


520.358 


4.86 


119 


46 


785.328 


8.46 


109 


45 


833.336 


4.86 


93 


52 



TABLE 7 

Absorption Lines in the Afterglow Spectrum of 
GRB 071003 



A 


z 


Transition 


W a 


a(W) b 


(A) 






(A) 


(A) 


3549.69 


0.37223 


Fell 2586 


< 2.51 




3568.06 


0.37223 


Fell 2600 


2.33 


0.59 


3837.72 


0.37223 


Mg II 2796 


2.48 


0.20 


3847.65 


0.37223 


Mg II 2803 


2.14 


0.19 


3915.45 


0.37223 


Mgl 2852 


1.02 


0.17 


4032.63 


1.60435 


CIV 1548 


0.22 


0.06 


4039.88 


1.60435 


CIV 1550 


< 0.28 




4351.92 


1.60435 


ALU 1670 


< 0.14 




5003.26 


1.10019 


Fell 2382 


0.20 


0.05 


5276.54 


1.60435 


Znll 2026 


< 0.08 




5399.79 


0.37223 


Call 3934 


0.61 


0.07 


5417.99 


0.93740 


Mg II 2796 


0.61 


0.05 


5432.79 c 


1.10019 


Fell 2586 


0.46 


0.05 


5447.85 


0.37223 


Call 3969 


0.46 


0.07 


5872.31 


1.10019 


Mg II 2796 


0.80 


0.05 


5888.27 


1.10019 


Mg II 2803 


0.68 


0.06 


6105.90 


1.60435 


Fell 2344 


< 0.17 




6206.91 


1.60435 


Fell 2382 


0.26 


0.04 


6240.46 


1.60435 


Fell* 2396a 


0.25 


0.04 


6265.95 


1.60435 


Fell* 2405 


< 0.16 




6282.68 


1.60435 


Fell* 2411b 


0.18 


0.03 


6284.57 






0.72 


0.12 


6734.47 






0.97 


0.15 


6737.28 


1.60435 


Fell 2586 


0.16 


0.04 


6772.60 


1.60435 


Fell 2600 


0.27 


0.05 


7301.58 


1.60435 


Mg II 2803 


0.17 


0.05 


7430.06 


1.60435 


Mgl 2852 


< 0.24 




8091.56 




0.92 


0.13 


8436. 10 d 






0.86 


0.26 


8534.91 d 






0.72 


0.17 


8599. 02 d 






1.34 


0.17 



a Equivalent widths are rest-frame values and assume 
the redshift given in Column 2. 
b Limits are 2a statistical values. 
c Blended with Mg II A2803 at z = 0.937. 
d These features may be residuals from sky subtrac- 
tion. 
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TABLE 8 

Optical Light-Curve Fits: Color Change 



Model Description 


Aj%-1 


A,0i (6 -a) 


A/3 -2 


fh 


x7" 


Fully monochromatic 
Uniformly chromatic bump 
Variably chromatic bump 
Chromatic rebrightening 
Chromatic bump+rebrightening 




0.22±0.27 
0.66±0.33 


0.75±0.33 





1.05±0.47 


1.09±0.47 






0.77±0.31 
0.84±0.31 


0.72±0.10 
0.68±0.10 
0.70±0.10 
1.26±0.11 
1.26±0.11 


125.765 / 81 
125.100 / 80 
120.339 / 79 
120.040 / 80 
113.713 / 78 



Note. — Summary of relevant parameters and y 2 f° r models allowing or disallowing color transi- 
tions and chromatic breaks between the various components. Values without uncertainties are fixed. 
Component is the fast-decay component, Component 1 is the bump, and Component 2 is the late 
rebrightening. The absolute late-time spectral index @2 is not a model parameter, but is fit externally 
after completion of the fit. 



TABLE 9 
Optical Light-Curve Fits: to 



Model Description 




dto 


dti 


dt 2 


X*/ v 






(s) 


(s) 


to 




Reference 













113.713 / 78 


Decay 




0.01 ± 3.01 








113.713 / 77 


Bump 







60.5±20.4 





112.700 / 77 


Bump (prompt pulse) 







125.0 





115.118 / 78 


Rebrightening 










1245±311 


111.149 / 77 



Note. — Summary of relevant parameters and \ 2 f° r models using a to 
different from the trigger time. In all cases, the favored color-change model 
(chromatic bump and rebrightening) was used. Values without uncertainties are 
fixed. Component is the fast-decay component, Component 1 is the bump, 
and Component 2 is the late rebrightening. 



TABLE 10 
Radio Modeling of GRB 071003 



.Parameter 


Value (broken power law) 


Value (unbroken) 


Value (unbroken w/ scintillation 11 ) 


Oib 


-0.11 ± 0.21 


0.27 ± 0.06 


0.34 ± 0.10 


Cta 


0.81 ± 0.25 






tbreak 


8.51 ± 3.78 








-1.11 ± 0.34 


-1.15 ± 0.44 


-1.15 ± 0.42 




2.45 / 4 


15.32 / 6 


6.07 / 6 



Note. — Best-fit parameters of a fit to the radio afterglow of GRB 071003 using a 
Bcucrmann ct al. (1999) broken power-law model versus an unbroken power-law model. The im- 
provement for the broken power-law fit is significant given the flux uncertainties, but due to inter- 
stellar scintillation may be coincidental. If a small amount of interstellar scintillation uncertainty 
is added in quadrature, an unbroken power-law fit is reasonable. 

a In this model, we added a 15% error to all X-band points and a 22% error to all C-band points. 



TABLE 11 
Model fluxes at i = 2.67 days 



Band/Filter 


E 


Flux 


U ncertainty 


cV 






X-ray 


1000 


0.036 


0.006 


u 


3.46 


3.17 


1.42 




2.55 


4.47 


0.30 


V 


2.25 


5.07 


0.27 


r 


1.97 


5.97 


1.14 


R 


1.88 


8.01 


0.47 


i 


1.61 


9.16 


0.56 


1 


1.54 


10.34 


0.54 


z 


1.34 


14.74 


2.29 


K' 


0.584 


33.59 


16.8 


X 


3.5e-5 


414.6 


91.8 


C 


2.0e-5 


256.1 


73.9 



Note. — Fluxes of the afterglow interpo- 
lated to t = 2.67 d after the BAT trigger us- 
ing all available X-ray, optical, and radio data. 
Galactic extinction (E(B — V) = 0.148 mag) 
is not accounted for; however, the X-ray flux 
is corrected for photoelectric absorption. 
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TABLE 12 

Extinction models for optical/X-ray fits of GRB 
071003 



model 




A 


V 






H v 











Y 2 Iv 


none 















_ 





.913 


± 





029 


12.4 / 8 


z = 0.372 


























Milky Way 


0. 


239 


± 





.093 


3 


.09 





.939 


± 





.028 


5.80 / 7 


SMC 





209 


± 





082 


2 


.74 


0. 


.934 


± 





.028 


5.95 / 7 


LMC 


0. 


.256 


± 





099 


3 


.41 





941 


± 





029 


5.87 / 7 


Calzetti 


0. 


279 


± 





.108 


4 


.05 





945 


± 





.029 


5.80 / 7 


z = 1.10 


























Milky Way 


0. 


133 


± 





.058 


3 


.09 





.935 


± 





.029 


7.16 / 7 


SMC 


0. 


127 


± 





.052 


2 


.74 





935 


± 





.028 


6.38 / 7 


LMC 





132 


± 





.057 


3 


.41 





934 


± 





.028 


7.16 / 7 


Calzetti 


0. 


247 


± 





.095 


4 


.05 





.957 


± 





.032 


5.78 / 7 


z = 1.60 


























Milky Way 





139 


± 





.048 


3 


.09 





.943 


± 





.028 


3.94 / 7 


SMC 





096 


± 





.037 


2 


.74 





934 


± 





.028 


5.77 / 7 


LMC 





131 


± 





.045 


3 


.41 





940 


± 





.028 


3.98 / 7 


Calzetti 


0. 


240 


± 





.093 


4 


.05 


0. 


.965 


± 





.033 


5.84 / 7 



Note. — Results of various fits to the contemporaneous 
optical and X-ray fluxes for extinction due to either the host 
galaxy or the intervening absorbers at z = 0.372 and z = 1.10. 
A small amount of extinction is required to accurately fit the 
data, but its nature is not strongly constrained. We adopt 
SMC-likc extinction at z = 0.372 in the discussion and plots 
based on the relative strength of the intervening absorber at 
this redshift in the spectrum. 



